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?  MISSILE2 ,  to  improve  its  accuracy  of  hinge-moment  prediction.  A 
method  based  on  strip  theory  and  shock-expansion  theory  is  used 
to  improve  the  prediction  of  fin  axial  center-of-pressure  posi¬ 
tion.  A  special  method  is  included  to  account  for  canard  and 
body  vortices  on  fin  axial  center-of-pressure  location.  Special 
extrapolation  procedures  for  the  underlying  data  base  were  added. 
Comparison  between  hinge-moment  data  and  prediction  show  good 
agreement  at  zero  roll  between  Mach  numbers  of  1.3  and  3.7,  but 
not  at  transonic  speeds.  When  body  or  canard  vortices  are  close 
to  tail  fins  on  the  leeward  side  of  the  body,  the  predictions 
were  usually  not  good.  Canard  fin  hinge  moments  were  generally 
well  predicted. 


PREFACE 


This  technical  report  covers  the  work  performed  under 
Contract  N00014-81-C-0267  from  March  1,  1981  to  November  30, 
1981.  The  program  is  sponsored  by  the  Office  of  Naval  Research 
Arlington,  Virginia.  Dr.  Robert  E.  Whitehead  is  the  Scientific 
Officer. 
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SUMMARY 


This  report  is  the  technical  report  on  the  first  phase  of  a 
two-part  investigation  of  the  hinge  moments  of  all-movable  con¬ 
trols  as  utilized  in  cruciform  missiles.  In  the  first  phase,  an 
attempt  has  been  made  to  develop  as  general  as  possible  a  hinge- 
moment  prediction  method  from  subsonic  to  hypersonic  speeds  for 
angles  of  attack  to  20°  and  control  deflections  up  to  15°.  The 
complete  range  of  roll  angles  has  been  included.  In  the  second 
phase,  a  study  is  to  be  carried  out  to  minimize  hinge  moments, 
and  an  experimental  investigation  is  to  be  undertaken  to  see  if 
the  low  hinge  moments  are  realized. 

Initially  a  literature  survey  was  made  to  collect  systematic 
data  on  hinge  moments  and  to  seek  methods  of  hinge-moment  predic¬ 
tion.  The  state  of  the  art  was  reviewed  with  respect  to  a  number 
of  factors  which  can  influence  hinge  moments  including  wing-alone 
nonlinearities,  airfoil  section  effects,  fin-fin  and  fin-body 
interferences,  gaps,  canard  and  body  vortices,  and  trailing-edge 
shock-wave,  boundary- layer  interaction.  One  effect  found  to  be 
of  significance  is  that  of  airfoil  thickness  ratio  and  thickness 
distribution  which  can  have  a  significant  effect  on  the  center  of 
pressure  of  fin  normal  force.  It  was  also  found  that  vortices 
such  as  those  of  the  body  or  canard  fins  passing  in  close  proximity 
to  the  all-movable  fin  can  significantly  change  its  center-of- 
pressure  position  and  normal  force  and  hence  its  hinge  moment. 
Sufficient  systematic  hinge-moment  data  were  found  to  use  for 
checking  the  results  of  any  prediction  method  for  the  effects  of 
angle  of  attack,  roll  angle,  Mach  number,  fin  planform,  and  fin 
airfoil  section(s).  An  unpublished  wing-alone  data  base  provided 
systematic  data  on  the  effects  of  fin  airfoil  section (s),  planform, 
angle  of  attack,  and  Mach  number  on  normal  force  and  axial  center- 
of-pressure  position. 

To  predict  hinge-moment  coefficient  requires  separate  pre¬ 
dictions  of  normal-force  coefficient  and  axial  center-of-pressure 
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location.  An  existing  missile  aerodynamic  force  and  moment  com¬ 
puter  program  (MISSILE  1  or  MISSILE  2)  generally  predicts  normal- 
force  coefficients  with  sufficient  accuracy,  but  not  axial 
center-of-pressure  locations.  An  equivalence  was  found  between 
the  axial  center  of  pressure  of  the  fin  on  body  with  that  of  the 
wing  alone  at  the  same  value  of  the  normal-force  coefficient. 

Thus  if  experimental  data  on  the  axial  center  of  pressure  of  the 
wing  alone  are  available,  they  can  be  used  in  the  hinge-moment 
prediction  method  provided  its  thickness  distribution  is  the  same 
as  that  of  the  fin.  For  those  cases  where  such  data  are  not 
available  but  where  the  leading-edge  shock  is  attached  viewed  in 
planes  normal  to  the  fin  planform  but  parallel  to  the  body  axis, 
a  method  was  found  for  predicting  axial  center-of-pressure  loca¬ 
tion.  The  method  uses  a  combination  of  strip  theory  and  shock- 
expansion  theory  in  a  separate  computer  program  to  obtain  axial 
center-of-pressure  location.  This  information  is  used  as  input 
to  MISSILE  2.  Also  a  method  for  determining  the  effects  of 
vortices  on  the  fin  hinge  moment  was  added  to  MISSILE  2.  With 
these  and  several  other  modifications  to  improve  accuracy,  a  ver¬ 
sion  of  MISSILE  2  to  predict  hinge-moment  coefficient  was  developed. 
This  program  is  applicable  to  fins  with  unswept  trailing  edges. 
However  a  means  for  treating  fins  with  swept  trailing  edges  is 
outlined. 

For  a  fin  of  general  planform  it  appears  possible  to  obtain 
accurate  hinge  moments  using  a  combination  of  the  program  DEMON  II 
and  the  strip  theory,  shock-expansion  program  referred  to  above. 
However,  no  example  is  carried  out. 

The  hinge-moment  version  of  MISSILE  2  has  been  applied  to  a 
variety  of  cases  to  determine  its  range  of  applicability  with 
reasonable  accuracy.  The  following  conclusions  are  borne  out  by 
comparison  of  experiment  and  theory. 

1.  The  method  has  been  successfully  applied  for  the  following 
range  of  Mach  number,  1.3  <  M^  <  3.7. 
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2.  While  the  method  appears  to  be  fairly  good  for  leeward 
canard  fins,  it  is  not  accurate  for  leeward  tail  fins  near  strong 
body  vortices  or  canard  vortices. 

3.  The  use  of  the  panel  normal  forces  obtained  from  DEMON  II 
can  increase  the  accuracy  of  hinge-moment  prediction  for  windward 
fins  if  used  rather  than  the  normal  forces  determined  from 
MISSILE  2. 

4.  For  zero  roll  angle  the  predicted  tail  fin  hinge  moments 
with  no  canard  deflections  are  good  and  those  with  canard  deflec¬ 
tions  are  fair.  At  a  roll  angle  of  45°,  the  predicted  tail  fin 
hinge  moments  are  not  good  with  or  without  canard  deflection 
because  the  canard  vortex  effects  on  the  tail  fin  are  not  well 
predicted  for  the  vortices  in  close  proximity. 

5.  The  method  is  not  accurate  at  transonic  speeds  principally 
because  no  method  exists  to  account  for  wing  section  effects  on 
axial  center-of-pressure  position. 
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1.  INTRODUCTION 


The  choice  of  fin  size  and  planform  shape  to  provide  a 
missile  with  adequate  control  but  small  hinge  moments  has 
been  a  difficult  problem  at  least  since  World  War  II.  The 
problem  has  generally  been  studied  by  means  of  expensive  and 
time-consuming  wind-tunnel  testing.  The  weights  of  the  con¬ 
trol  actuators  and  their  power  supplies  are  substantial. 

While  some  attention  has  been  given  to  analytical  means  for 
predicting  hinge-moments  for  missiles,  no  attempt  seems  to 
have  been  made  to  develop  general  hinge-moment  prediction 
technology  and  to  apply  it  to  the  search  for  low  hinge-moment 
all-movable  controls  (fins) .  It  is  the  purpose  of  this  report 
to  describe  a  preliminary  method  for  estimating  hinge  moments 
of  all-movable  controls.  The  Mach  number  range  of  interest 
is  0.5  <  M^  <  4.0. 

Many  missiles  use  all-movable  controls  for  trim  and 
maneuver  because  the  greatest  control  effectiveness  is  usually 
obtained  by  utilizing  all  the  potentially  available  fin  area 
for  control.  Preliminary  design  methods  for  determining  the 
characteristics  of  all-movable  controls  have  been  available 
only  for  the  linear  range  of  angle  of  attack  and  control 
deflection.  As  missiles  have  been  pushed  to  higher  angles  of 
attack  in  search  of  greater  maneuverability,  the  required 
control  deflections  have  also  increased.  The  nonlinearities 
resulting  from  higher  angle-of-attack  operation  have  made  the 
prediction  of  the  control  characteristics  even  more  difficult. 
Accordingly,  the  normal  forces  and  hinge  moments  developed 
by  such  controls  have  usually  been  measured  in  wind  tunnels 
during  the  design  stage  to  assure  that  they  will  produce  the 
desired  trim  and  maneuverability  within  the  hinge-moment 
capabilities  of  the  control  actuators. 
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Recent  studies  (ref.  1)  of  the  forces  acting  on  the  fins 
of  cruciform  wing-body  combinations  at  high  angles  of  attack 
have  had  some  success  in  predicting  the  high  angle  of  attack 
nonlinearities.  There  is  therefore  now  some  prospect  for 
developing  a  preliminary  method  for  estimating  hinge-moments 
of  all-movable  controls  at  high  angles  of  attack.  It  is  the 
purpose  of  the  present  report  to  present  such  a  method.  The 
method  is  preliminary  because  further  work  is  planned  under 
the  contract,  and  the  method  is  constructed  so  that  it  can 
utilize  a  better  data  base  which  is  planned  for  the  future. 

This  work  is  being  conducted  under  Contract  N00014-81-C- 
0267  from  the  Office  of  Naval  Research  which  is  supported  by 
the  Flight  Dynamics  Laboratory  of  the  Air  Force  Wright 
Aeronautical  Laboratories.  The  work  is  being  carried  out  in 
two  phases  as  follows. 

Phase  I  -  Development  of  an  analytical  capability  for  pre¬ 
dicting  hinge  moment  for  aerodynamically-tailored  fins,  and 
validation  and  refinement  of  the  analytical  method  following 
correlation  of  the  results  from  the  predictive  method (s) 
with  existing  missile  fin  hinge-moment  data. 

Phase  II  -  An  unconstrained  look  at  approaches  for  lower¬ 
ing  fin  hinge  moment  while  retaining  missile  control 
effectiveness  shall  be  made  using  the  analytical  capability 
developed  during  Phase  I.  Upon  analysis  of  the  results,  not 
less  than  two  and  not  more  than  six  planforms  reflecting  the 
most  promising  approaches  shall  be  selected  and  tested.  The 
experimental  data  shall  be  analyzed  for  the  purpose  of  refining 
the  analytical  capability  further  to  improve  the  accuracy  in 
predicting  hinge  moments  of  aerodynamically  tailored,  low 
hinge-moment  fins. 

In  this  report  we  start  with  a  general  discussion  of  the 
problem  followed  by  a  summary  of  the  state  of  the  art.  Then 
the  predictive  method  is  described.  A  series  of  comparisons 
between  experiment  and  theory  are  carried  out  and  discussed. 


In  the  concluding  remarks,  the  successes  and  shortcomings  of 
the  method  will  be  covered. 

2.  GENERAL  CONSIDERATIONS 

2.1  General  Observations 

The  hinge  moment  depends  on  the  normal  force  developed 
by  a  fin  together  with  the  center-of-pressure  position  of 
the  normal  force  with  respect  to  the  hinge  line.  A  number 
of  methods  exist  for  determining  the  fin  normal  force  for 
various  ranges  of  the  relevant  parameters.  However, 
sufficiently  accurate  methods  for  determining  longitudinal 
center-of-pressure  position  do  not  generally  exist.  This 
void  is  addressed  herein. 

In  the  estimation  of  hinge-moment  coefficient  and  its 
comparison  with  data,  care  must  be  taken  in  assessing  the 
accuracy  of  prediction.  For  a  control  with  small  hinge 
moments,  the  percent  difference  between  the  theory  and  data  is 
not  a  good  measure  of  accuracy  since  if  the  hinge-moment  were 
zero,  the  error  would  be  percentagewise  infinite  even  for 
excellent  agreement.  A  better  measure  is  how  much  the  differ¬ 
ences  between  experiment  and  theory  represents  in  terms  of  the 
hinge-moment  lever  arm  expressed  as  percent  of  the  root  chord. 
Based  on  the  present  state  of  the  art,  if  the  error  represents 
less  than  about  2  percent  of  the  root  chord,  it  is  considered 
that  good  agreement  is  obtained. 

2.2  Parameters  of  Interest 

Some  discussion  of  the  ranges  of  the  parameters  to  be 
considered  will  help  fix  ideas.  First  by  an  all-movable  control 
we  mean  that  the  entire  fin  rotates  about  some  lateral  axis 
which  projects  from  the  missile  body.  The  emphasis  will  be  on 
the  cruciform  case  rather  than  the  planar  case,  although  at  <f>  =  0 
they  are  essentially  the  same.  Also  for  small  fin-fin  interference 


they  do  not  differ  substantially.  Ratios  of  body  diameter  to 
wing  span  from  0  to  well  above  0.5  are  covered. 

With  regard  to  fin  parameters,  we  based  the  aspect  ratio 
on  two  opposing  fins  joined  together  at  their  root  chords. 
Aspect  ratios  from  1/2  to  2.0  are  covered  in  the  data  base 
used  in  the  method  for  the  normal  force  of  the  fin  mounted  on 
the  body.  Extrapolation  up  to  aspect  ratio  approaching  3  is 
possible.  A  future  data  base  up  to  3*  =  4  and  =  4.5  is 
planned.  Taper  ratios  of  0  to  1  are  included.  For  the  most 
part  we  will  be  dealing  with  fins  with  unswept  trailing  edges 
in  the  data  base.  However,  a  means  of  specifying  an  equivalent 
fin  with  an  unswept  trailing  edge  will  be  given  for  fins  with 
swept  trailing  edges.  The  method  will  be  structured 
sufficiently  broadly  that  fins  of  polygonal  planform  with 
many  straight  edges  can  be  handled.  For  all-movable  controls 
the  aspect  ratio  range  of  importance  is  generally  greater 
than  unity. 

Besides  the  geometric  parameters  of  the  basic  configura¬ 
tion,  we  are  interested  in  the  following  ranges  of  the  param¬ 
eters  : 

0°  <  a  <  20°;  angle  of  attack 
-15°  <  6  <  15°;  control  deflection  angle 
0.5  <  <  4.5;  Mach  number 

0°  <  <t>  <  180°;  roll  angle 

Most  of  our  effort  will  be  devoted  to  all-movable  controls 
not  operating  in  the  wake  of  another  set  of  cruciform  surfaces. 
The  computer  program  will  be  applicable  to  the  case  of  all¬ 
movable  controls  operating  in  the  wake  of  forward  fins,  but 
the  nonlinearities  encountered  are  more  numerous  than  for 
forward  controls  because  of  the  wake  of  the  forward  controls 
and  that  generated  by  the  afterbody  behind  the  forward  surfaces 
We  will  make  some  comparisons  between  experiment  and  theory  for 
this  case  to  determine  our  ability  to  handle  it. 


2.3  Factors  Influencing  Hinge  Moment 

A  number  of  factors  have  an  influence  on  the  hinge-moment 
coefficients  of  all-movable  controls.  A  list  is  given  below 
as  a  guide  to  the  discussion  of  the  state  of  the  art  in  the 
following  section. 

1.  Nonlinearities  in  the  wing-alone  normal-force  curve 
(versus  a) 

2.  Nonlinearities  in  the  wing-alone  longitudinal  center- 
of-pressure  curve  (versus  a) 

3.  Fin-body  interference  (no  body  vortices) 

4.  Fin-body  interference  due  to  body  vortices 

5.  Fin-fin  interference 

6.  Canard  vortices  (effects  on  aft  surfaces) 

7.  Adverse  fin-body  interference 

8.  Fin  trailing-edge  shock-wave /boundary- layer  inter¬ 
action 

9.  Gap  effects 

10.  Aeroelastic  effects 

11.  Fin  choking 

3.  STATE  OF  THE  ART 
3.1  Introductory  Remarks 

A  literature  survey  has  been  made  through  Defense 
Technical  Information  Center  (DTIC)  and  NASA  Scientific  and 
Technical  Information  Facility  to  find  references  regarding 
all-movable  controls.  The  references  found  have  been  classified 
under  the  following  headings. 


A.  Gap  Effects 

B.  Work  at  Ames  Research  Center 

C.  Work  at  Langley  Research  Center 

D.  General 

The  bibliography  is  presented  in  Appendix  A.  The  overall 
methods  for  determining  the  characteristics  of  all-movable 
controls  are  generally  fluid  dynamic  computer  programs  which 
determine  the  entire  flow  field  with  the  exception  of 
reference  2  which  presents  an  engineering  method.  The  present 
method  is  a  considerable  extension  of  that  method. 


3.2  Nonlinearities  in  the  Wing-Alone 
Normal-Force  Curve 

One  of  the  inputs  to  the  calculation  method  is  the  wing- 
alone  normal- force  curve.  There  are  numerous  theories  for 
determining  wing  normal-force  coefficient,  but  they  are  not 
adequate  to  determine  the  wing  normal  force  to  high  angles  of 
attack  with  confidence  for  general  planforms.  It  is  not  our 
intent  to  try  to  remedy  this  general  defect  of  theory,  but 
instead  we  will  rely  on  data  base  methods. 

Previously,  the  only  wing-alone  data  base  for  the  normal- 
force  coefficients  from  low  to  high  Mach  numbers  and  to  high 
angle  of  attack  of  wings  which  systematically  vary  in  aspect 
ratio  and  taper  ratio  was  obtained  on  a  reflection  plane  and 
contained  unknown  effects  of  the  plate  boundary  layer.  Wing- 
alone  normal-force  curves  have  been  "constructed"  from  existing 
data  for  aspect  ratios  of  0.5,  1.0,  2.0;  taper  ratios  of  0, 

0.5,  and  1.0;  and  Mach  numbers  of  0.8,  1.2,  2.0,  and  3.0  for 
use  in  reference  1. 

Recently  Stallings  and  Lamb,  reference  3,  published  a 
wing-alone  data  base  for  Mach  numbers  from  1.6  to  4.50  based 
on  pressure  distribution  models.  The  set  of  wings  tested  by 
Stallings  and  Lamb  have  the  geometric  characteristics  shown  in 
Figure  1.  It  is  noted  that  the  aspect  ratio  goes  up  to  4. 


These  airfoil  sections  are  not  of  uniform  thickness  ratio 
but  are  of  constant  thickness  because  of  the  necessity  of 
housing  the  pressure  tubing.  The  values  of  normal-force 
coefficient  agree  closely  enough  with  those  of  reference  1 
for  the  range  tested  so  that  no  changes  are  required  on 
reference  1. 

Recently  the  Stallings-Lamb  wings  were  tested  in  the  Ames 
6-  by  6-Foot  Supersonic  Wind  Tunnel  at  M^  =  0.6,  0.8,  0.9,  and 
1.2  to  complete  the  Mach  number  range.  The  data  have  been 
reduced  and  some  of  them  are  utilized  in  the  present  report. 
This  work  was  accomplished  under  Contract  DAAG29-79-C-0020 
for  the  Army  Research  Office,  Durham,  N.C.  with  support  from 
ARO,  NAVAIR,  and  MICOM.  The  models  for  the  supersonic  testing 
were  designed  and  constructed  by  the  Langley  Research  Center 
and  lent  to  Nielsen  Engineering  &  Research,  Inc.  (NEAR)  for 
testing  in  the  Ames  Research  Center  6-  by  6-Foot  Wind  Tunnel. 
Dr.  Robert  E.  Singleton,  Director,  Engineering  Sciences 
Division,  AROD,  has  given  us  permission  to  use  some  of  these 
unpublished  data  in  this  report. 

A  wing-alone  data  base  due  to  Baker  (ref.  4)  was  obtained 
on  a  reflection  plane  using  semispan  models.  The  aspect- 
ratio  range  was  1/2  to  2  and  the  Mach  number  range  up  to 
M^  =  3.  At  high  angles  of  attack  and  Mach  numbers,  the  data 
are  known  to  exhibit  loss  of  normal  force  of  unknown  extent 
due  to  boundary-layer  effects  of  the  reflection  plane.  This 
data  base  will  be  termed  the  Fidler-Baker  data  base  hereafter 
since  the  wings  tested  by  Baker  are  similar  to  wings  tested 
earlier  by  Fidler. 

3.3  Nonlinearities  in  the  Wing-Alone  Longitudinal 
Center-of-Pressure  Position 

The  longitudinal  center-of-pressure  position  determines 
the  lever  arm  for  the  hinge  moment.  Any  variation  of  this 
position  with  angle  of  attack  for  the  wing  alone  would  be 


expected  to  also  be  present  to  some  degree  for  the  fin  in  the 
presence  of  the  body. 

The  theoretical  determination  of  the  wing-alone  center- 
of-pressure  location  is  in  the  same  state  of  development  as 
the  prediction  of  the  wing-alone  normal  force.  We  must  rely 
on  data  to  a  large  extent.  With  this  in  mind  we  examine  the 
question  of  the  available  systematic  data  on  the  matter. 

However,  first  consider  what  symmetrical  thickness  can  do  to 
the  center-of-pressure  location  of  a  two-dimensional  super¬ 
sonic  airfoil. 

It  is  known  that  a  symmetric  thickness  distribution 
causes  the  airfoil  center-of-pressure  location  at  supersonic 
speed  to  move  forward  from  the  zero  thickness  position  by  an 
amount  that  increases  with  thickness  ratio  and  Mach  number. 

In  fact,  it  is  possible  to  use  Busemann  second-order  theory 
or  shock-expansion  theory  to  show  this. 

Using  both  of  these  theories,  calculations  have  been  made 
for  the  effect  of  thickness  on  center-of-pressure  positions  for 
the  three  airfoils  of  Figure  2.  These  airfoils  have  the  same 
thickness  ratio  but  vary  in  the  length  of  the  region  of  maxi¬ 
mum  thickness.  Without  any  thickness,  the  center-of-pressure 
position  denoted  by  x/c  is  at  the  midcord  of  the  airfoil,  and 
with  thickness  is  moved  forward  an  amount  (Ax/c)t.  The  values 
of  this  latter  quantity  are  given  in  Figure  3  for  the  three 
airfoils  as  calculated  by  the  two  theories  for  =  4.5.  It 
is  noted  that  large  forward  shifts  occur  which  must  be  accounted 
for  if  x/c  is  to  be  predicted  to  within  say  0.02.  Note  that  the 
Busemann  theory  and  shock-expansion  theory  agree  at  a  =  0°,  but 
shock-expansion  theory  exhibits  an  angle-of-attack  effect. 

Also  the  double-wedge  airfoil  shows  the  least  effect  of  thick¬ 
ness.  Figure  4  shows  the  effects  of  Mach  number  for  one  of 
the  airfoils,  and  demonstrates  how  (Ax/c)t  increases  with  M^. 
These  effects  occur  when  the  leading-edge  shock  wave  is  attached 


The  foregoing  results  show  why  center-of-pressure  data 
from  a  set  of  wings  of  systematic  planforms  do  not  yield  system¬ 
atic  center-of-pressure  locations  unless  the  airfoil  sections 
are  uniform  or  vary  systematically  wing  to  wing.  They  also 
show  that  a  systematic  hinge-moment  data  base  requires  system¬ 
atic  airfoil  sections. 

There  exist  two  sets  of  wings  varying  systematically  in 
aspect  ratio  and  taper  ratio  for  which  longitudinal  center-of- 
pressure  data  have  been  obtained  to  large  angles  of  attack; 
the  set  shown  in  Figure  1  and  the  Fidler-Baker  set  shown  in 
Figure  5.  It  is  noted  that  the  two  sets  of  wings  have  dis¬ 
similar  airfoil  sections,  and  the  airfoil  section  changes  from 
wing  to  wing  within  each  set.  A  comparison  of  thickness  ratios 


is  mad^  in 

the  following 

tabulation 

between  the  two 

sets 

of 

wings  for 

wings  of  the  same  aspect 

ratio  and  taper 

ratio. 

Stallings-Lamb 

Fidler-Baker 

Wing 
(fig.  1) 

(t/cr)root 

Wing 
(fig.  5) 

(t/cr) root 

m 

A 

P2 

.044 

T31 

.025 

0.5 

0.5 

P4 

.042 

T14 

.025 

1.0 

0 

P5 

.062 

T15 

.037 

1.0 

0.5 

P6 

.083 

Til 

.037 

1.0 

1.0 

P8 

.088 

T23 

.049 

2.0 

0.5 

The  Fidler-Baker  wings  were  tested  on  a  boundary-layer  reflec¬ 
tion  plate,  and  induced  plate  boundary-layer  separation  at 
large  values  of  a  which  invalidated  the  data.  The  Stallings- 
Lamb  wings  were  tested  as  pressure-distribution  wings,  a  fact 
accounting  for  their  greater  thickness  ratios.  Comparisons 
between  x/cr  for  the  two  sets  of  wings  will  now  be  made  for 
=  0.8,  1.2,  2.0,  and  3.0. 

In  Figures  6-13  comparisons  will  be  made  between  the 
normal-force  coefficients  and  axial  center-of-pressure  locations 
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for  the  five  wing  planform  common  to  the  Fidler-Baker  and 
Stallings-Lamb  data  bases  for  Mach  numbers  of  0.8,  1.2,  2.0, 
and  3.0.  Before  the  presentation  of  these  comparisons  some 
points  with  regard  to  the  data  are  noteworthy.  First  the 
data  for  the  wings  of  the  Fidler-Baker  data  base  were  cor¬ 
rected  for  the  tares  in  normal- force  and  hinge-moment 
coefficients  at  a  =  0  before  determining  axial  center-of- 
pressure  position,  x/cr>  The  Fidler-Baker  wings  were  not 

tested  at  M  =1.2  but  were  tested  at  M  =1.15  and  1.3.  The 
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results  presented  for  =  1.2  represent  a  linear  interpola¬ 
tion  between  the  data  for  these  two  Mach  numbers.  No  data 
for  the  T14  wing  exist  for  Mk  =  0.8  or  1.2,  but  the  P4  data 
are  included  for  completeness. 

Consider  now  the  comparison  of  the  wing  normal-force  curves 
at  M^  =  0.8  given  in  Figure  6.  The  first  apparent  result  is 
that  the  two  sets  of  wings  have  quite  different  normal- force 
curves  at  =  0.8  as  a  result  of  different  wing  sections. 

There  is  a  difference  in  their  stall  characteristics,  the 
thinner  wings  of  the  Fidler-Baker  data  base  generally  stalling 
at  higher  angles  of  attack  for  aspect  ratios  of  1  or  more. 

Consider  now  the  values  of  x/cr  for  the  same  wings  in 
Figure  7.  The  wings  of  the  Fidler-Baker  data  base  show  large 
rearward  travel  of  the  center-of-pressure  as  angle  of  attack 
increases  in  all  cases.  The  wings  P6  and  P8  do  not  exhibit 
such  large  rearward  shifts.  The  two  sets  of  wings  show 
generally  much  different  x/cr  characteristics. 

Figures  6  and  7  stow  the  large  effects  of  wing  section 
on  CN  and  x/cr  which  have  implications  with  regards  to  pre¬ 
dicting  hinge  moment.  Since  a  general  data  base  including 
planform  and  wing-section  effects  is  not  available,  we  must 
confine  our  predictive  capaoility  to  wings  with  existing  data 
or  must  rely  on  theory.  However,  a  theory  which  accounts  for 
wing-section  effects  on  x/cr  does  not  exist  for  the  angle  of 
attack  in  which  we  are  interested.  Possibly  transonic  full 


potential  theory  would  work  at  the  lower  angles,  but  its 
application  is  expensive  and  limited.  A  general  prediction 
method  for  hinge  moment  around  =  0.8  thus  requires  more 
basic  work. 

We  will  now  consider  the  state  of  affairs  for  M  =1.2. 

00 

Figure  8  presents  the  comparisons  for  the  normal-force  coef¬ 
ficients.  Generally  the  effects  of  wing  section  on  normal- 
force  coefficients  are  small  up  to  wing  stall.  The  thinner 
wings  of  the  Fidler-Baker  data  base  exhibit  no  significant 
stall  while  P 5,  Pg,  and  Pg  do. 

The  x/cr  comparisons  of  the  wings  for  =  1.2  are  given 
in  Figure  9.  The  thinner  wings  generally  have  centers  of 
pressures  near  the  DATCOM  values  for  small  angles  of  attack, 
which  move  rapidly  backward  with  increasing  angle  of  attack. 

For  the  Pg,  Pg,  and  Pg  wings  which  exhibited  normal-force  stall, 
the  centers  of  pressure  are  in  front  of  those  for  the  corre¬ 
sponding  Fidler-Baker  wings  above  the  stall.  At  low  angles  of 
attack  the  center-of-pressure  positions  are  aft  of  those  of 
the  Fidler-Baker  wings. 

We  have  a  problem  at  Mro  =  1.2  with  regard  to  predicting 
center-of-pressure  positions  at  high  angle  of  attack  because 
of  the  lack  of  a  suitable  theory.  For  low  angles,  where  the 
leading-edge  shock  may  be  attached,  we  will  subsequently 
present  a  method. 

Consider  now  the  comparison  of  the  normal-force  curves  for 
the  two  sets  of  wings  at  Moo  =  2.0  as  given  in  Figure  10.  The  curves 
are  in  very  good  agreement  except  at  the  highest  angles  of  attack 
where  wings  T31  and  T14  exhibit  stall  at  angles  less  than  50°. 

The  corresponding  center-of-pressure  locations  are  given  in 
Figure  11.  The  points  in  the  stalled  region  for  T14  and  T31  have 
been  marked  with  flags  since  the  stall  may  be  induced  by  reflec¬ 
tion  plate  boundary  layer.  It  is  noted  that  the  Stallings-Lamb 
wings  generally  exhibit  center-of-pressure  locations  in  front 
of  those  for  the  Fidler-Baker  wings,  a  result  that  would  be 


expected  in  the  range  of  shock  attachment  because  of  their 
greater  thickness.  This  difference  would  also  be  expected  near 
a  *  0°,  but  this  is  not  the  case  for  T31  and  T14.  This  result 
may  be  due  to  inaccuracies  in  determining  x/cr  near  a  =  0°. 

For  the  7R  = 1,  A  =  1  wings  larger  differences  are  shown  at 
small  a  than  at  large  a.  This,  again,  may  be  an  accuracy 
problem. 

The  normal-force  curves  and  the  axial  center-of-pressure 
curves  for  the  two  sets  of  wings  at  M  =3.0  are  shown  in 
Figures  12  and  13,  respectively.  The  Fidler-Baker  wings  T31, 
T14,  T15,  and  Til  show  earlier  stall  than  the  corresponding 
Stallings-Lamb  wings.  This  phenomenon  is  reflected  in  the 
axial  center-of-pressure  positions  shown  in  Figure  13  with 
flagged  symbols  in  the  neighborhood  of  a  =  50°  and  above. 

These  flagged  data  points  are  suspect.  The  observations 
made  concerning  the  differences  in  the  results  for  the  two 
sets  of  wings  at  =  2.0  also  hold  generally  at  =  3.0. 

The  fairly  regular  differences  shown  give  hope  of  predicting 
these  differences. 

3.4  Fin-Body  Interference;  No  Body  Vortices 

3.4.1  Introductory  remarks.-  A  fin  in  the  presence  of 
the  body  can  act  differently  than  it  does  when  part  of 
the  wing  alone  because  of  the  nonuniform  flow  field  induced 
at  the  fin  location  by  the  body  and  because  the  round  body 
does  not  present  an  infinite  reflection  plane.  We  now  con¬ 
sider  the  effects  of  fin-body  interference  without  the 
complicating  effects  of  body  vortices  which  are  taken  up  in 
the  next  section. 


We  are  interested  in  the  effects  of  fin-body  interference 
on  the  quantities  CnW(B)  and  x^^  since  both  of  these 
quantities  are  used  in  calculating  hinge  moment.  The  fin-body 
interference  depends  on  the  following  parameters  in  a 
complicated  way. 


a 


angle  of  attack 
Mm;  free-stream  Mach  number 
<p ;  roll  angle 

a/sm;  radius/semispan  ratio 

n;  number  of  fins  (equally  spaced) 

6^,62/63,64;  fin  deflection  angles 
JR;  fin  aspect  ratio 
X;  fin  taper  ratio 

We  will  confine  our  attention  to  n  =  4. 

Consider  now  what  theories  are  available  for  determining 
CnW(B)  and  (*/cr ^ w (b) ' 

3.4.2  Theory.-  For  subsonic  speeds  several  theories  are 
available  for  determining  the  fin  load  distribution.  First 
there  is  the  slender-body  theory  of  Spreiter  and  Sacks 
(ref.  7)  for  6  =  0,  Adams  and  Dugan  (ref.  8)  for  6^0.  There 
are  also  various  lifting-line  theories  and  lifting-surface 
theories.  These  methods,  however,  are  generally  worked  out 
for  zero  wing  thickness  and  apply  to  the  linear  range  of 
angle  of  attack  characteristics.  A  special  vortex-lattice 
program  has  been  developed  by  Dillenius  (ref.  9)  which  applies 
to  empennages  of  quite  general  geometric  configuration.  This 
program,  applicable  to  Mtt  =  0,  treats  multiple  fins  which  are 
not  alike  and  which  can  be  unequally  spaced  around  the  fuselage 

For  supersonic  speeds,  there  are  a  number  of  linear 
methods  that  can  be  used  including  quasi-cylindrical  theory 
(ref.  10)  and  paneling  methods.  Of  this  latter  type,  the 
PANAIR  program  (ref.  11)  is  an  example  as  well  as  DEMON  II 
(ref.  12) ,  which  is  a  panel  program  especially  adapted  to 
cruciform  missiles.  For  supersonic  speeds,  DEMON  II  can  be 
used  for  fins  of  arbitrary  planforms.  Codes  based  on  the 
Euler  equations  have  been  used  for  supersonic  wing-body  com¬ 
binations  for  supersonic  axial  Mach  numbers.  The  code  of 


Soloman,  Wardlaw  and  Baltakis  (ref.  13)  and  that  of  Klopfer 
and  Nielsen  (ref.  14)  are  examples. 

Generally  the  panel  codes  applicable  to  wing-body 
combinations  require  significant  amounts  of  computer  time. 

They  are  applied  to  unusual  planforms  or  for  calculation  of 
benchmark  cases.  For  preliminary  design  we  will  rely  as  much 
as  possible  on  a  data  base  or  on  engineering  methods  for 
determining  both  CNw(b)  and  (x/cr)w(B). 

3.4.3  Experiment.-  A  data  base  for  CnW(B)  and  (X/cr)w^Bj 
should  cover  a  range  of  wing  planforms,  angles  of  attack,  roll 
angles,  Mach  numbers,  and  deflection  angles.  Such  a  data  base 
does  not  exist  for  variable  6.  For  6=0  the  Fidler-Baker 
data  base  exists  from  which  only  the  CNW(B)  results  have  been 
incorporated  into  Program  MISSILE  (Refs.  1  and  2) .  The  data 
base  for  (x/cr)w^Bj  is  not  systematic  since  the  airfoil 
sections  varied  among  wings  of  the  set. 

In  reference  2  it  was  shown  that  (X/cr)w^  was  closely 
approximated  by  (x/cr)w  for  the  wing  alone.  The  correspondence 
of  the  center  of  pressure  of  the  wing  alone  and  the  fin  in  the 
presence  of  the  body  is  based  on  an  equivalence  assumption. 

One  way  of  putting  the  assumption  is  that  at  the  same  value  of 
normal-force  coefficient  based  on  planform,  the  longitudinal 
and  lateral  center-of-pressure  positions  of  the  wing  alone  and 
the  fin  are  the  same.  In  other  words,  the  fin  is  assumed  to 
exhibit  the  same  nonlinearities  as  the  wing  at  a  given  value 
of  fin  normal-force  coefficient. 

We  will  now  compare  the  axial  center-of-pressure  positions 
of  the  wings  alone,  (x/cr)w,  as  a  function  of  fin  normal-force 
coefficient  (ref.  5)  with  those  for  the  fins  in  the  presence 
of  the  body,  (x/cr)w^B^,  using  data  (ref.  15)  for  the  wings  of 
the  Fidler-Baker  data  base.  The  first  such  wing  results  are 
shown  in  Figure  14  for  Mach  numbers  of  0.8,  1.2,  2.0,  and 
3.0  for  AR  =0.5  and  X  =  0.5.  Some  explanation  of  the  results 
are  necessary.  The  data  for  the  fin  in  the  presence  of  the 


body,  taken  from  reference  15,  are  given  for  various  roll 
positions  on  the  windward  side  of  the  body  to  minimize  any 
effects  of  body  vortex  interference.  They  have  not  been 
corrected  for  zero  shifts  of  normal  force  or  hinge  moment, 
and  therefore  contain  some  degree  of  uncertainty.  The  wing- 
alone  data  have  been  taken  from  reference  5,  and  they  have 
been  corrected  for  zero  shifts. 

For  all  Mach  numbers  it  is  noted  that  the  data  for  x/cr 
for  the  fin  at  different  values  of  <}>  tend  to  correlate  in  a 
small  band.  Also  the  wing-alone  data  generally  goes  through 
this  band  within  about  ±0.02.  The  nonlinear  trend  of  the 
center  of  pressure  with  angle  of  attack  measured  for  the  wing 
alone  is  thus  duplicated  by  the  fin  in  the  presence  of  the 
body.  It  is  also  apparent  that  the  equivalence  principle 
described  above  does  a  good  job  of  correlating  the  data  for 
different  roll  angles  and  angles  of  attack. 

The  results  for  four  more  wings  are  given  in  Figures  15 
to  18  in  accordance  with  the  following  table: 


Configuration 
(ref.  15) 

_A_ 

Figu 

B1T14 

1 

0 

15 

B1T15 

1 

1/2 

16 

BlTll 

1 

1 

17 

B1T23 

2 

1/2 

18 

In  Figure  15  for  B1T14,  wing-alone  data  are  not  available  for 
T14  at  M  =0.8  and  1.2  so  that  only  the  data  for  the  fin  in 
the  presence  of  the  body  are  shown  for  this  combination*.  A 
mean  line  has  been  drawn  through  the  data  for  these  two  Mach 

_ 

While  a  wing  of  comparable  planform  was  tested  at  these  Mach 
numbers  in  the  Stallings-Lamb  data  base,  it  has  a  different 
thickness  distribution  from  that  for  T14. 


numbers  to  simulate  the  missing  wing-alone  results.  It  is 
noted  for  this  delta  fin  that  the  center-of-pressure  varia¬ 
tions  with  angle  of  attack  and  Mach  number  are  not  large, 
possibly  making  it  a  candidate  for  a  low  hinge  moment  all¬ 
movable  control. 

Figure  16  for  B1T15  exhibits  an  interesting  phenomenon 
at  =  0.8.  For  a  range  of  CN  between  0.87  and  0.97,  x/cr 
of  the  wing  alone  is  triple  valued  for  a  given  CN<  It  is  not 
clear  from  the  data  for  the  fin  in  the  presence  of  the  body 
whether  this  nonlinearity  is  present  for  the  fin.  We  have 
compared  the  data  for  the  wing  alone  with  that  of  the  fin  at 
two  different  roll  angles  in  Figure  19.  Note  that  at  the 
horizontal  position,  the  fin  duplicates  the  wing  alone 
nonlinearity  past  stall.  For  <j>  =  30°  the  fin-on-body  dupli¬ 
cates  the  wing  alone  up  to  the  stall.  To  predict  fin  center- 
of-pressure  beyond  the  stall  from  wing  alone  data  in  the 
transonic  range  is  thus  not  promising. 

Figure  17  for  B1T11  shows  a  slight  stall  for  the  /R  =  1, 

X  =  1  fin  at  Mw  =  0.8  but  it  appears  not  to  be  important.  The 
large  center-of-pressure  travel  for  7R  =1,  X  =  1.0  and  that 
for  /R  =  2,  X  =  0.5  shown  in  Figures  17  and  18  at  =  0.8  and 
1.2  make  it  difficult  to  obtain  small  hinge-moment  coefficients 
with  them  in  the  transonic  range. 

It  appears  that  the  wing-alone  axial  center-of-pressure 
position  is  thus  a  good  approximation  to  that  for  the  fin  in 
the  presence  of  the  body  for  the  five  wings  of  the  Fidler- 
Baker  data  base  for  the  Mach  number  range  =  0.8  to  3.0. 
However,  in  cases  of  severe  wing-alone  stall,  the  approxima¬ 
tion  breaks  down. 

3.5  Fin-Body  Interference  Due  to  Body  Vortices 

3.5.1  Introductory  remarks.-  A  body  of  revolution  above 
an  angle  of  attack  depending  principally  on  its  length-to- 
diameter  ratio  will  generate  a  pair  of  symmetric  vortices  on 


its  leeward  side.  At  a  still  larger  angle  of  attack  the 
symmetric  pair  can  become  asymmetric  particularly  for  subsonic 
flow.  The  present  method  of  predicting  hinge  moments  could 
be  applied  to  the  asymmetric  vortex  case  if  we  knew  the 
vortex  positions  and  strengths.  A  rough  rule  of  thumb  is  that 
asymmetric  vortices  do  not  develop  if  the  crossflow  Mach 
number  is  greater  than  about  0.5  (ref.  16).  This  fact  taken 
together  with  the  fact  that  asymmetric  vortices  form  at  angles 
of  attack  above  about  20°-30°  for  body  of  moderate  length-to- 
diameter  ratio  leads  to  the  plot  shown  as  Figure  20.  The 
present  method  has  not  been  applied  to  the  region  marked 
"Subcritical ,  asymmetric." 

In  the  region  above  the  curve  Mc  =  Mm  sin  a  =  0.5,  there 
are  generally  shock  waves  in  the  crossflow  plane  which  appear 
to  emanate  from  the  body  surface.  These  shock  waves  tend  to 
eliminate  the  asymmetric  vortices.  The  vortical  regions  are 
not  so  concentrated  as  for  subcritical  crossflow  Mach  numbers, 
but  tend  to  be  elongated  elliptical  regions. 

Body  vortices  are  generally  characterized  for  engineering 
purposes  by  their  strengths  and  their  positions.  Their  strengths 
are  given  by  a  nondimensionalized  value  of  the  circulation 
about  a  contour  including  the  vorticity.  The  vortex  position 
is  characterized  by  the  location  of  the  core  of  the  vortex  in 
the  crossflow  plane.  This  position  can  usually  be  deduced 
from  flow  surveys,  vapor- screen  photographs,  and  sometimes 
schlieren  photographs.  A  description  of  the  vortices  in  terms 
of  strengths  and  positions  is  the  first  step  in  determining 
the  effect  of  the  body  vortices  on  the  fin.  The  second  step 
is  determining  what  effect  these  vortices  have  on  fin  forces 
and  moments  using  aerodynamic  theory. 

3.5.2  Theory.-  We  consider  first  the  available  theories 
for  determining  the  vortex  strengths  and  positions.  These 
theories  fall  into  two  classes,  vortex  tracking  methods  and 


Euler  equations  methods.  Vortex  tracking  methods  consider 
the  crossflow  planes  of  a  body  of  revolution  and  introduce 
point  vortices  at  the  body  separation  line,  as  calculated  or 
as  input  empirically.  The  calculation  proceeds  downstream 
with  the  shed  vortices  following  the  calculated  streamlines. 

The  method  of  Mendenhall,  et  al.  (refs.  17-19)  is  an  example 
of  this  type  of  theory.  A  second  method  makes  use  of  the 
Euler  equations  and  known  positions  of  the  body  primary  vortex 
shedding  lines.  By  the  application  of  appropriate  boundary 
conditions  at  the  separation  line,  the  entire  flow  field  is 
calculated  in  a  supersonic  marching  procedure.  The  method 
of  Nielsen  and  Klopfer  (ref.  20)  is  the  only  present  example 
of  this  approach. 

We  next  consider  methods  for  calculating  the  effect  of 
the  vortices  on  the  fin  aerodynamics.  These  methods  include 
reverse-flow  methods  (ref.  1),  panel  methods  (ref.  12),  and 
Euler  methods  (ref.  20).  In  the  reverse-flow  method,  the 
upwash  distribution  induced  at  the  fin  by  the  external 
vortices  and  their  images  is  first  determined.  Then  a  span- 
wise  integral  of  this  upwash  distribution  multiplied  by  a 
weighting  factor  is  used  to  determine  the  fin  normal  force 
due  to  the  vortices.  Another  weighting  factor  yields  the 
lateral  center  of  pressure.  This  method  is  used  in  Program 
MISSILE  (ref.  1) .  However,  the  method  as  originally  formulated 
is  not  suitable  for  determining  hinge  moments,  and  an  engineer¬ 
ing  means  has  to  be  developed  where  none  now  exists. 

Panel  methods  are  applicable  to  determining  hinge  moments 
due  to  vortices.  The  use  of  Program  DEMON  II  (ref.  12)  in 
this  connection  will  be  demonstrated  subsequently.  The  Euler 
equation  method  of  determining  the  effect  of  body  vortices  on 
the  fin  loading  was  recently  developed  by  NEAR  and  is  a  finite- 
difference  method  not  presently  suitable  for  preliminary  design. 

3.5.3  Experiment.-  Consider  fist  the  experimental  data 
with  respect  to  strengths  and  positions  of  symmetric  body 


vortices  above  bodies  of  revolution  at  angles  of  attack. 

Vapor- screen  photographs  and  conventional  pitch- yaw  probes 
were  used  in  the  1950's  to  obtain  vortex  positions,  and  the 
vortex  strengths  were  inferred  principally  from  normal  forces 
or  normal-force  distributions.  These  data  have  been  correlated 
for  use  in  preliminary  design  by  several  investigators.  A 
correlation  of  these  data  by  Mendenhall  and  Nielsen  given  in 
reference  17  has  been  used  in  Program  MISSILE. 

In  recent  years  additional  measurements  in  this  area 
have  been  made  by  Oberkampf  (ref.  21)  using  conventional 
probes.  Schwind,  et  al.  (refs.  22  and  23)  and  Yanta,  et  al. 
(ref.  24)  have  also  made  measurement  using  laser-doppler 
velocimetry . 

With  regard  to  measuring  the  effect  of  vortices  on  a 
wing  pressure  distribution,  two  studies  have  been  made  at  low 
speed  by  Smith  and  Lazeroni  (ref.  25)  and  McMillan,  et  al. 

(ref.  26).  The  latter  authors  found  that  reverse-flow  theory 
provides  loads  of  sufficient  accuracy  compatible  with 
economy.  To  obtain  better  accuracy  probably  requires  account¬ 
ing  for  the  effect  of  the  fin  on  the  vortex  path  in  a  fully 
mutual  interference  calculation. 

3.6  Fin-Fin  Interference 

3.6.1  Introductory  remarks.-  Fin-fin  interference  occurs 
when  the  aerodynamic  characteristics  of  a  given  fin  are  changed 
as  the  result  of  the  presence  of  another  fin  of  the  set.  There 
is  also  fin-fin  interference  if  the  aerodynamic  characteristics 
of  a  given  fin  are  influenced  by  control  deflection  of  any  other 
member  of  the  set.  It  is  clear  that  fin-fin  interference  can 
influence  all-movable  control  effectiveness  and  hinge  moment 
and  must  be  accounted  for  in  a  hinge-moment  prediction  method. 

A  severe  case  of  fin-fin  interference  can  result  in  serious 
control  cross-coupling.  In  this  case  a  pitch  control  command 
might  result  in  yaw  control  and/or  roll  control  as  well.  It 


is  clear  that  fin-fin  interference  is  more  severe  for  cruci¬ 
form  fins  than  for  planar  fins  simply  because  of  the  closer 
proximity  of  the  neighboring  fins. 

Fin-fin  interference  occurs  when  one  fin  falls  within  the 
zone  of  influence  of  another  fin  as  at  subsonic  speed.  At 
supersonic  or  hypersonic  speed,  it  is  possible  that  the  fins 
are  not  within  each  others  fields  of  influence  so  that  no 
fin-fin  interference  occurs. 

3.6.2  Theory.-  Several  methods  exist  for  calculating 
fin-fin  interference  due  to  control  deflection  and  the  simplest 
of  these  is  probably  reverse-flow  theory.  Appendices  A  and  B 
of  reference  1  developed  this  theory  using  a  combination  of 
slender-body  theory  and  reverse-flow  theorems.  The  method 
determines  the  effect  of  control  deflection  on  fin  normal 
force  and  lateral  center  of  pressure,  but  is  not  adaptable 
to  determining  the  effect  on  longitudinal  center-of-pressure 
position.  Fin-fin  interference  for  no  control  deflection 
has  already  been  accounted  for  in  the  data  base  of  Program 
MISSILE. 

Panel  methods  such  as  Program  DEMON  (ref.  12) ,  specially 
constructed  for  cruciform  missiles,  can  be  used  to  study  fin- 
fin  interference.  Euler  codes  can  be  developed  to  handle 
general  fin-fin  interference  in  cruciform  missiles  for  the 
case  of  supersonic  marching. 

Theoretical  fin-fin  interference  factors  for  unit  fin 
deflections  have  been  calculated  using  Program  MISSILE  and  are 
presented  in  Table  1  of  reference  1.  These  factors  are 
applicable  to  fins  which  are  wholly  within  each  others  zones 
of  influence.  For  a  fin  partially  within  the  field  of 
influence  of  other  fins,  the  percentage  of  the  fin  planform 
within  the  field  of  influence  is  determined.  This  percentage 
is  applied  to  the  above  mentioned  fin-fin  interference  factors 
as  a  multiplication  factor.  The  method  gives  the  right 


limit  when  a  fin  is  not  within  the  field  of  influence  of  a 
deflected  fin. 

3.7  Canard  Vortices;  Effects  on  Aft  Surfaces 

The  influence  of  canard  vortices  on  the  hinge  moments  of 
all-movable  tail  controls  is  complicated.  For  angles  of  attack 
which  are  relatively  small  and  for  zero  roll  angle,  the  hori¬ 
zontal  canard  fins  produce  vortices  which  may  be  sufficiently 
close  to  the  tail  fins  to  influence  their  hinge  moments.  At 
high  angles  of  attack  (<j>  =  0°)  the  canard  vortices  can  pass 
high  enough  above  the  tail  fins  to  have  a  small  effect  on 
their  hinge  moments.  However,  at  high  angles  of  attack  the 
afterbody  between  the  canards  and  the  empennage  can  shed  its 
own  vorticity  which  interacts  with  the  tail  fins  and  has  a 
strong  effect  on  their  hinge  moments.  The  interaction  between 
the  canard  vortices  and  the  afterbody  will  influence  the 
strengths  and  positions  of  both  the  afterbody  and  canard 
vortices. 

The  foregoing  effects  encountered  at  zero  roll  angle  are 
made  more  complicated  at  nonzero  roll  angle  by  vortices  from 
all  four  canard  panels  which  interact  with  each  other  and  with 
the  afterbody  in  an  asymmetric  way.  To  account  for  these 
vortex  effects  on  tail  hinge  moments  requires  a  knowledge  of 
the  vortex  strengths  and  positions  at  the  tail  axial  location 
and  then  an  application  of  a  method  for  determining  the  tail 
fin  reactions.  The  strengths  and  positions  of  the  canard 
vortices  at  the  canard  trailing  edge  can  be  obtained  from 
some  knowledge  of  the  span  loading  distribution;  then  the 
canard  vortices  must  be  tracked  over  the  afterbody,  the  origin 
of  the  afterbody  vortices  determined,  and  the  strengths  and 
positions  of  the  afterbody  vortices  determined  accounting  for 
their  interaction  with  the  canard  vortices. 

MISSILE  2  (ref.  2)  contains  a  vortex-tracking  routine 
to  take  account  of  the  above  phenomena.  In  principle  it  can 


be  applied  to  the  problem  of  determining  the  vortex  strengths 
and  positions  at  the  tail  location.  The  method  subsequently 
to  be  described  for  determining  the  effects  of  vortices  on  fin 
hinge  moments  can  then  be  used  to  determine  the  tail  fin  hinge 
moments.  To  be  able  to  account  accurately  for  all  nonlinearities 
on  tail  fin  hinge  moments  is  probably  more  than  can  be  expected 
for  the  present  method,  but  some  comparisons  will  be  made  to 
determine  the  state  of  the  art. 


3.8  Adverse  Fin-Body  Interference 


At  supersonic  speeds  at  moderate  angles  of  attack  and  at 
hypersonic  speeds  at  low  angles  of  attack  significant  nonlinear 
effects  occur  on  the  fin  normal  force  due  to  the  interference 
of  the  body.  Account  of  these  nonlinearities  must  be  taken 
in  determining  hinge  moments. 


The  fin  in  the  presence  of  the  body  at  6  =  0°  and  <p  =  0° 
has  a  normal-force  coefficient  Cnw^  based  on  its  planform 
area,  and  the  wing  alone  at  the  same  angle  of  attack  has  a 
normal  force  coefficient  Cnw«  The  interference  of  the  body 
on  the  fin  is  measured  by  a  factor  Kw  given  by 


KW  = 


CnW(B) 


'N, 


w 


This  factor  is  normally  greater  than  unity  at  low  speeds  and 
low  angles  of  attack  because  of  body  upwash,  denoting  favorable 
interference.  The  experimental  values  of  these  factors  for 
an  aspect  ratio  1  delta  wing  on  a  body  are  shown  in  Figure  21. 
The  body  vortex  effects  have  been  taken  out  of  Cnw(b)  by 
theoretical  means.  it  is  seen  that  the  favorable  interference 
at  low  Mach  numbers  and  low  angles  of  attack  can  become  very 
unfavorable  at  high  Mach  numbers  and  high  angles  of  attack. 


The  favorable  interference  manifest  at  low  Mach  numbers 
and  angles  of  attack  is  known  to  be  the  result  of  body  upwash 
which  increases  the  fin  normal  force.  However,  at  high  Mach 


numbers  and  angles  of  attack,  the  net  interference  is  a  result 
of  changes  in  flow  upwash  angle,  dynamic  pressure,  and  Mach 
number  at  the  positions  to  be  occupied  by  the  fin  in  the 
body  flow  field.  The  effect  of  compressibility  and  the  bow 
wave  on  these  quantities  basically  controls  the  value  of  Kw. 

A  fair  estimate  of  the  fin  normal  force  can  be  made  by 
strip  theory  if  the  flow  field  in  the  position  to  be  occupied 
by  the  fin  is  known.  Euler  codes  exist  which  will  calculate 
this  flow  field  quickly  for  supersonic  marching.  However, 
since  Program  MISSILE  uses  a  data  base  for  cN^(b)'  the  values 
of  Kw  are  known  experimentally  from  the  data  base.  They  are 
used  in  extending  the  data  base  to  configurations  with  a/sm 
ratios  different  from  0.5,  that  of  the  data  base. 

3.9  Boundary-Layer,  Shock-Wave  Interaction 

At  an  unswept  (or  supersonic)  sharp  trailing  edge  of  a 
supersonic  wing  on  the  upper  surface,  the  flow  is  above  free- 
stream  speed,  and  must  traverse  an  oblique  shock  to  come  back 
to  free-stream  static  pressure  behind  the  wing.  Under  certain 
circumstances  the  shock  moves  forward  from  the  trailing  edge 
and  separates  the  boundary  layer  approaching  the  trailing 
edge.  The  results  presented  by  Ferri,  reference  27,  for  a 
biconvex  airfoil  with  semiapex  angle  of  llo20'  tested  at 
M^  =  2.13  at  Re  =  710,000  are  interesting.  They  show  effects 
of  boundary- layer  separation  at  a  =  0°  over  the  last  15  percent 
of  the  chord  and  at  a  =  20°  over  the  last  50  percent  of  the 
chord. 

If  boundary- layer ,  shock-wave  interaction  occurs  on  the 
upper  surface,  the  pressure  coefficients  on  the  top  rear  part 
of  the  wing  are  increased,  so  that  the  normal  force  per  unit 
area  is  decreased.  As  a  result  the  center-of-pressure  loca¬ 
tion  should  move  forward  and  change  the  hinge  moment. 

The  pressure  distributions  taken  by  Stallings  and  Lamb 
on  their  wings  sheds  some  light  on  the  extent  of  separated  flow 


on  the  upper  surface  of  the  wings  near  the  trailing  edge. 

For  this  purpose  the  pressure  distributions  on  the  JR  =  2 
rectangular  wing  have  been  examined.  This  wing  has  15°  semi¬ 
wedge  angles  at  the  leading,  side,  and  trailing  edges.  It  was 
tested  at  a  Reynolds  number  of  slightly  over  700,000  based  on 
the  wing  chord.  The  orifices  numbered  52,  57,  and  63  in 
Figure  22  have  the  following  distances  from  the  trailing 
edge : 


Orifice 

x'  /cr 

52 

0.177 

57 

0.118 

63 

0.059 

The  measured  pressure  distributions  on  the  upper  and  lower 
surfaces  of  the  wing  at  a  =  10°  are  compared  in  Figure  23(a) 
with  the  predictions  of  shock-expansion  theory.  Note  that  on 
the  upper  surface  the  last  three  orifices  show  uniform  pressure 
on  the  beveled  trailing  edge.  The  pressure  is  slightly  higher 
than  predicted  by  shock-expansion  theory.  Some  separation 
from  the  shoulder  to  the  trailing  edge  seems  to  be  present. 

At  any  rate  the  difference  between  experiment  and  theory 
accounts  for  about  0.015  in  x/cr.  The  data  in  Figure  23(b) 
for  a  =  20°  probably  exhibit  separation.  No  comparison  is 
made  with  theory  since  the  angle  of  attack  is  above  that  for 
shock  detachment. 

Data  and  theory  are  compared  in  Figures  23(c)  and  (d) 

for  M  =4.60  for  a  =  20°.  There  is  little  or  no  evidence 

00 

of  trailing-edge  separation.  These  data  show  that  the 
pressure  coefficient  on  the  upper  surface  of  a  wing  at  high 
angle  of  attack  and  Mach  number  is  so  close  to  the  free-stream 
value  that  boundary- layer ,  shock-wave  interactions  on  the  upper 
surface  approaching  the  trailing  edge  could  not  have  a  signifi¬ 
cant  effect  on  normal-force  or  center-of-pressure  position. 
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We  will  generally  ignore  the  effects  of  trailing-edge, 
shock-wave,  boundary-layer  interaction  on  center  of  pressure 
in  the  hinge-moment  predictive  method.  This  should  prove  to 
be  a  fair  approximation  for  wedge-shaped  trailing  edges.  How¬ 
ever,  on  the  basis  of  the  results  of  Ferri  the  same  assumption 
should  not  be  made  for  wings  with  airfoil  sections  curved  at 
the  trailing  edge  such  as  biconvex  sections. 

3.10  Gap  Effects 

3.10.1  Introductory  remarks.-  A  gap  exists  between  an 
all-movable  control  and  the  body  on  which  it  is  mounted.  Two 
kinds  of  gaps  occur  in  practice  which  are  worth  differentiating 
For  a  fin  at  6  =  0  there  is  generally  a  streamwise  gap  of  nearl; 
uniform  width  between  the  fin  and  the  body.  This  is  the  kind 
of  gap  which  is  generally  modeled  in  theory.  In  practice 

the  fin  is  rotated  about  a  hinge  line  so  that  a  gap  of  non- 
uniform  width  occurs  between  the  body  and  the  fin  for  6  #  0. 

For  very  large  deflections,  the  fin  can  even  project  above  the 
body.  To  treat  this  kind  of  nonuniform  gap  analytically  is 
clearly  very  difficult  but  it  is  amenable  to  panel  methods  with 
some  complexity.  However,  for  small  control  deflections  the 
gap  is  fairly  well  approximated  by  a  uniform  gap. 

Principally  as  a  result  of  theoretical  studies  the 
effect  of  gaps  on  the  fin  characteristics  have  at  times  thought 
to  be  large.  On  the  basis  of  potential  theory,  the  span 
loading  at  the  gap  falls  to  zero  from  the  substantial  value 
it  has  for  no  gap.  However,  it  appears  that  for  gap  sizes 
used  in  practice,  much  of  the  adverse  effect  predicted 
theoretically  is  overcome  by  viscosity.  That  is  not  to  say 
that  for  nonuniform  gaps  due  to  6  there  are  not  substantial 
effects.  Let  us  examine  the  theoretical  basis  for  gap  effects 
as  well  as  the  experimental  results. 

3.10.2  Theory.-  There  are  two  theoretical  methods  that 
have  been  pursued  in  the  study  of  gap  effects.  The  first 


approach  is  that  of  Bleviss  and  Struble  (refs.  28  and  29)  who 
use  linearized  supersonic  wing  theory  to  study  the  effect  of 
a  gap  between  a  wing  panel  and  a  vertical  reflection  plane. 

The  second  approach  based  on  slender-body  theory  has  been  used 
by  Mirels  (ref.  30)  and  Dugan  and  Hikido  (ref.  31)  to  study 
slender  wings  mounted  on  bodies  of  revolution  with  nonuniform 
gaps.  Bleviss  and  Struble  find  that  the  presence  of  a  gap, 
as  small  as  0.05  of  the  exposed  semispan,  can  cause  losses  in 
lift  as  large  as  40  percent  of  that  for  no  gap  (ref.  28,  Fig.  4) 
Dugan  and  Hikido  also  find  large  losses  of  lift  even  for 
minute  gaps.  Those  theoretical  results  are  considerably 
modified  by  viscosity. 

3.10.3  Experiment.-  Several  experimental  investigations 
of  gap  effects  are  available  (refs.  32  and  33) .  Careful 
systematic  measurements  were  reported  by  Dugan  (ref.  32)  on 
the  effects  of  a  gap  on  the  lift  developed  by  an  aspect  ratio 
2.31  delta  wing  mounted  on  a  body  of  revolution  with  a  body 
diameter  to  span  ratio  of  0.216  for  zero  gap.  The  tests  were 
carried  out  at  a  Mach  number  of  1.4  and  a  Reynolds  number  of 
slightly  less  than  4  x  10^  based  on  the  wing  root  chord.  Two 
cases  were  tested,  the  effect  of  a  at  6  =  0  and  the  effect  of 
6  at  a  =  0.  Some  results  are  now  quoted  from  Figure  5  of  the 
report. 

For  a  gap/wing  semispan  ratio  of  0.025,  the  ratio  of  lift 
with  a  gap  to  that  with  no  gap  was  0.75  at  a  =  1°,  0.95  at 
a  =  5°,  and  0.97  at  a  =  10°.  For  a  =  0  but  variable  <$ ,  the 
comparable  results  are  0.72  at  6  =  1°,  0.94  at  6  =  5°,  0.97 
at  a  =  10°.  At  very  small  angles  of  attack  or  control  deflec¬ 
tion  significant  loss  of  lift  occurred;  by  the  time  a  reaches 
5°  and  6  reaches  10°  very  little  loss  occurs. 

An  investigation  (ref.  34)  of  an  aspect  ratio  4  all¬ 
movable  delta  wing  on  a  body  of  revolution  also  contains  some 
data  on  gap  effects.  The  results  quoted  here  are  for 
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gap/wing  semispan  ratios  of  0.0057  and  0.023,  a  Mach  number  of 
1.46,  and  a  Reynolds  number  of  5.7  *  10^  based  on  a  body  length 
of  10.5  inches  with  a  root  chord  of  2.25  inches.  It  was  con¬ 
cluded  that  over  the  range  of  gap/semispan  ratios  tested 
there  was  little  or  no  effect  on  lift,  drag,  or  hinge  moment 
at  M  =  1.46. 
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It  thus  appears  that  for  gap/semispan  ratios  up  to  about 
0.025  significant  loss  of  fin  lift  occurs  only  at  very 
small  values  of  a  or  6 .  At  higher  values  of  a  or  6,  the  loss 
disappears  so  that  there  is  generally  no  basic  loss  in  maneuver 
capability.  If  larger  gaps  are  used,  larger  losses  occur 
(ref.  33) .  These  results  are  not  confirmed  for  high  supersonic 
or  hypersonic  Mach  numbers. 

3.11  Aeroelastic  Effects,  Fin  Choking 

Several  other  factors  can  influence  hinge  moments,  and 
they  are  mentioned  here  for  the  sake  of  completeness.  The 
first  factor  to  be  considered  is  aeroelasticity .  If  fins  are 
made  as  thin  as  possible  to  reduce  weight  and  are  still  subject 
to  large  loads  in  order  to  produce  high  maneuverability,  then 
significant  deformation  of  fins  under  load  is  to  be  expected. 

A  study  to  determine  to  what  extent  the  normal  forces  and 
hinge  moments  of  an  all-movable  control  are  changed  as  a  result 
is  needed.  In  particular,  it  appears  that  some  reduction  in 
hinge  moments  could  be  made  possible  by  aeroelastic  tailoring. 
The  determination  of  aeroelastic  effects  on  hinge  moments 
will  not  be  taken  into  account  in  the  present  prediction  method. 

The  second  factor  which  can  influence  hinge  moments  is 
"fin  choking"  which  can  occur  under  special  circumstances.  In 
this  phenomenon  the  axial  flow  between  the  two  lower  fins  of 
a  cruciform  wing-body  combination  in  the  x  roll  position 
suddenly  chokes  by  means  of  a  shock  normal  to  the  axial  direc¬ 
tion.  A  high  pressure  region  exists  between  the  fins  and  the 
lines  of  constant  pressure  are  nearly  lateral.  An  example 
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of  the  pressure  field  in  this  case  is  given  in  Figure  25  of 
reference  35  based  on  data  provided  by  Mr.  R.  L.  Stallings 
and  Mr.  W.  J.  Monta  of  Langley  Research  Center,  NASA.  The 
conditions  under  which  this  example  of  fin  choking  occurred 
are  <)>  =  45°,  a  =  50°,  and  Ma  =  2.7.  The  limits  of  operation 
to  avoid  choking  have  not  been  determinined. 

4.  PREDICTIVE  METHOD  FOR  ESTIMATING  HINGE  MOMENTS 
4.1  Introductory  Remarks 

A  method  for  estimating  hinge  moments  for  the  range  of 
angle  of  attack,  0°-20°,  and  the  range  of  control  deflection 
±15°,  is  not  generally  available.  A  method  applicable  to  the 
linear  range  is  given  in  reference  36.  Without  going  into  the 
details  of  that  method,  it  can  be  said  that  the  present  method 
is  an  extension  of  that  method  to  increase  its  range  of  appli¬ 
cability.  It  is  our  desire  that  the  present  method  include 
effects  of  roll  angle  and  be  applicable  from  subsonic*  to  high 
supersonic  speeds.  It  will  also  turn  out  that  the  method  we 
are  adopting  will  include  within  its  framework  the  possibility 
of  determining  the  hinge  moments  for  all-movable  tail  controls 
affected  by  wing-tail  interference.  It  is  also  an  objective 
to  treat  wing  planforms  of  quite  general  shape  so  long  as  these 
can  be  represented  by  connected  straight  line  segments. 

Calculating  hinge  moment  involves  determining  the  fin 
normal-force  coefficient  and  its  longitudinal  center-of- 
pressure  position.  We  have  available  an  engineering  tool  for 
predicting  fin  normal  force  coefficients  (MISSILE  2,  ref.  2) 
but  it  does  not  predict  axial  center-of-pressure  location 
with  sufficient  accuracy  for  hinge-moment  purposes.  Accord¬ 
ingly,  our  approach  is  to  use  MISSILE  2  as  the  starting  point 
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It  will  turn  out  that  the  present  method  needs  further  work 
for  application  at  transonic  speeds. 


of  the  prediction  technique  and  to  make  the  necessary  changes 
to  it  to  improve  hinge-moment  prediction. 


There  are  a  number  of  changes  or  additions  required  to 
MISSILE  2,  and  a  list  follows. 


Provide  accurate  wing-alone  center-nf-cressure  posi¬ 
tions  as  a  function  of  angle  of  attack  as  input. 


Add  the  capability  for  determining  the  effect  of 
free  vortices  on  fin  axial  center-of-pressure 


position. 


Add  a  better  means  for  extrapolation  outside  the 
/R  = 2  limit  of  the  data  base  of  MISSILE  2. 


Add  a  better  means  of  extrapolation  above  M  =  3  in 


the  data  base. 


Change  the  method  of  accounting  for  fin-fin  inter¬ 


ference. 


Determine  a  method  for  handling  fins  with  swept 
trailing  edges  within  the  framework  of  MISSILE  2, 


Item  1  comes  about  because  the  most  accurate  available 


means  for  obtaining  (x/cr)w  is  needed  to  predict  hinge 
moments.  The  most  accurate  means  available  differs  from  case 
to  case  so  the  calculation  is  done  external  to  the  main  program 
Item  2  is  needed  because  no  engineering  method  for  carrying 
out  this  task  has  hitherto  existed.  Items  3  and  4  are  needed 
because  much  of  the  available  data  for  evaluating  the  method 
lie  above  A R  =  2  and  M^  =  3.0.  Also  much  of  the  region  of 
interest  lies  above  these  values.  Item  5  changing  the  method 
of  accounting  for  fin-fin  interference  is  incorporated  into 
the  program  for  technical  reasons  which  will  be  more  important 
in  the  future.  Item  6  is  included  to  increase  the  configura¬ 
tional  range  of  applicability  of  the  program. 
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In  the  following  sections,  we  will  describe  these  above 
items.  We  will  not,  however,  give  a  full  description  of 
MISSILE  2  since  this  information  is  covered  in  references  1 
and  2 . 

4.2  Wing-Alone  Center-of-Pressure  Locations 

4.2.1  Introductory  remarks.-  The  wing-alone  longitudinal 
center-of-pressure  location  is  not  solely  a  function  of  angle 
of  attack,  Mach  number,  and  planform,  but  depends  also  on  the 
airfoil  section  or  sections  used  by  the  wing.  There  are 
several  options  in  determining  the  curve  of  (x/cr)w  versus  a 
for  the  wing  alone,  and  this  information  is  generated  external 
to  the  computer  program  for  convenience.  The  first  and  best 
option  is  to  have  experimental  data  for  the  wing  alone  (with 
the  correct  wing  thickness  distribution)  and  use  this  directly 
In  case  such  data  are  not  available,  we  have  several  options 
open  to  us.  If  the  leading-edge  shock  wave  is  attached  at 
any  angle  of  attack  when  considering  planes  normal  to  the  wing 
planform  and  parallel  to  the  body  axis,  then  a  combination  of 
linear  theory  (DATCOM) ,  shock-expansion  theory  and  data  is 
used  to  establish  the  curve.  In  case  the  shock  wave  is  not 
attached  at  any  angle  of  attack,  only  special  cases  can  be 
handled.  Let  us  consider  the  attached  shock  case  first. 

4.2.2  Attached  shock  case.-  A  brief  description  of  the 
method  is  given  before  a  detailed  explanation  of  how  it  is 
carried  out.  First,  the  center  of  pressure  for  an  infinitely 
thin  wing  is  determined  from  DATCOM  or  elsewhere,  and  a  cor¬ 
rection  for  thickness  is  obtained  using  shock-expansion  theory 
This  establishes  the  center-of-pressure  location  up  to  the 
shock  detachment  angle  of  attack,  ct(jet.  For  a  =  45°  it  is 
assumed  that  wing  section  has  little  effect  on  x/cr,  and  an 
empirical  correlation  of  this  quantity  is  used  to  establish 
its  value  at  a  =  45°.  Between  a.  .  and  a  =  45°,  the  curve  is 
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faired  by  hand  using  the  data  for  the  closest  wing  of  the 
data  base  as  a  guide.  We  will  be  considering  only  wings  with 
symmetric  thickness  distributions. 


To  determine  if  the  leading-edge  shock  of  the  wing  is 
attached,  we  need  only  consider  the  a  =  0  case.  Let  the 
nose  semi-apex  angle  in  the  streamwise  direction  be  0n,  and 
let  the  semi-vertex  angle  of  a  wedge  for  shock  detachment 
be  0det.  Then  if 


We  have  chosen  to  apply  shock-expansion  theory  to  planes 
parallel  to  the  body  axis  since  the  flow  near  the  fin-body 
juncture  is  principally  in  this  direction.  For  a  swept-back 
wing  of  high  aspect  ratio,  simple  sweep  theory  dictates  use  of 
simple  sweep  theory  in  planes  normal  to  the  wing  leading  edge. 
However,  such  wings  are  not  usually  used  for  all-movable 
controls. 

Consider  now  the  determination  of  the  wing  longitudinal 
center-of-pressure  position  for  the  shock-attached  case.  The 
first  step  is  the  determination  of  the  center-of-pressure 
position  for  the  wing  with  no  thickness  using  linear  theory. 
For  most  planforms  this  quantity  is  conveniently  obtained  from 
DATCOM.  For  delta  wings  the  center  of  pressure  is  always  at 
the  2/3  root  chord.  For  rectangular  wings  with  effective 
aspect  ratios,  8®  ,  of  unity  or  greater 
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For  wings,  the  leading  and  trailing  edges  of  which  are  composed 
of  straight  line  segments,  a  panel  method  such  as  DEMON  II, 
reference  12,  can  be  used.  If  time  and  cost  permits,  DEMON  II 
can  be  used  to  determine  the  fin  center-of-pressure  in  the 
presence  of  the  body,  which  can  be  used  rather  than  that  of 
the  wing  alone  for  greater  accuracy  of  prediction. 

The  second  step  in  determining  x/cr  is  to  correct  the 
lifting-surface  theory  value  for  wing  thickness  effects.  This 
correction  gives  a  forward  movement  of  the  center-of-pressure 
position  which  varies  with  angle  of  attack. 

Consider  now  the  change  in  center-of-pressure  location 
due  to  wing  thickness  distribution.  Let  the  wing  have  a 
straight*  trailing  edge,  and  consider  the  following  sketch. 


A  given  streamwise  section  as  shown  will  be  considered  to  act 
at  the  given  angle  of  attack,  a,  and  the  given  Mach  number  M^. 
By  shock-expansion  theory  the  section  will  have  a  longitudinal 
center-of-pressure  location  x/c  where  x  is  measured  from  the 
leading  edge  and  c  is  the  local  chord.  It  will  also  have 
section  lift  coefficient  c^(y).  For  no  thickness,  the  local 
section  center  of  pressure  will  be  at  x/c  =  0.5. 

Accordingly,  the  contribution  to  the  moment  due  to 
thickness  is  proportional  to  cc^(x  -  c/2).  The  shift  in  x 
for  the  entire  wing  is  using  strip  theory 
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To  carry  out  the  calculations  based  on  Equation  (5) ,  a 
computer  program  has  been  constructed.  First  the  section 
normal-force  coefficient  and  section  center-of-pressure  posi¬ 
tions  are  calculated  at  specified  spanwise  locations  on  the 
wing.  The  equations  to  carry  out  these  calculations  for  the 
airfoil  sections  by  shock-expansion  are  given  in  Appendix  B. 
It  is  assumed  that  the  airfoil  sections  are  symmetrical  and 
that  the  surface  can  be  approximated  by  a  connected  series 
of  straight  line  segments.  The  computer  program  then 
carries  out  the  calculation  indicated  by  Equation  (5) . 
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At  this  point  we  compare  the  measured  center-of-pressure 
locations  of  Stallings-Lamb  wings  (ref.  3)  with  the  predictions 
of  the  above  method  for  the  angle  of  attack  and  Mach  number 
ranges  of  shock  attachment.  Comparisons  are  made  for  10  wings 

from  M  =  1.60  to  4.60  as  follows: 
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The  degree  of  correlation  between  experiment  and  prediction 
is  now  systematically  reviewed. 

Considering  first  Figure  24,  it  is  noted  that  three  sets 
of  information  for  each  Mach  number  are  given;  the  data,  the 
linear  theory,  and  the  linear  theory  plus  thickness  correction. 
The  angle  of  attack  for  shock  detachment  increases  as  the  Mach 
number  increases.  The  thickness  correction  increases  as  Mach 
number  increases.  Considering  Figures  24-26  for  all  three 
delta  wings,  the  prediction  with  thickness  correction  and  the 
data  are  in  satisfactory  agreement.  Forward  shifts  of  the 
center-of-pressure  position  up  to  6  percent  of  the  root  chord 
are  exhibited. 

Consider  now  the  results  for  A  =  0.5  wings  in  Figures  27 
through  30.  For  JR  =0.5  the  agreement  is  not  good  at  all. 

DATCOM  predicts  positions  generally  in  front  of  the  experi¬ 
mental  center-of-pressure  position,  and  the  thickness  correction 
further  aggravates  the  agreement.  It  is  clear  that  some 


phenomenon  not  accounted  for  by  the  theory  is  involved  for 

the  JR  =0.5,  A  =  0.5  wing  which  did  not  act  for  the  JR  =0.5, 

A  =  0  wing.  There  is  a  strong  component  of  nonlinear  normal 
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force  proportional  to  a  due  to  leading-edge  and  side-edge 
separation  which  could  act  well  back  on  the  wing  and  cause  the 
observed  effect.  This  possiblity  was  not  pursued  in  an 
attempt  to  verify  it  or  to  improve  the  prediction  method  for 
several  reasons.  First,  JR  =  0.5  controls  are  not  common. 
Secondly,  for  a  wing  of  the  same  planform  but  of  different 
thickness  distribution,  we  can  use  the  data  of  Figure  27  and 
correct  it  for  thickness  distribution.  The  data  of  Figure  28 
for  the  JR  =1.0,  A  =  0.5  wing  exhibit  the  same  disagreement 
between  experiment  and  theory  for  M^  =  1.6  and  2.5,  but  at 
the  higher  Mach  number  the  agreement  is  much  better.  The 
JR  =2,  A  =  0.5,  and  JR  =4,  A  =  0.5  wings  shows  good  agreement 
everywhere  in  the  attached  shock  regime. 

Consider  now  the  A  =  1.0  wings  in  Figures  31-33.  The 
JR  =0.5  wing  shows  some  disagreement  at  Mro  =  2.86  but  is  good 
at  the  higher  Mach  numbers.  The  JR  = 1  and  2  wings  show  good 
agreement. 

In  the  discussion  we  have  considered  two  approaches: 

Method  A:  DATCOM  plus  thickness  corrections 

Method  B:  Use  the  data  with  corrections  to  account 
for  any  change  in  thickness  distribution. 

For  wings  of  various  taper  ratio,  the  method  is  now  summarized. 
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Above  the  angles  of  shock  detachment  from  the  wing  leading 
edge,  the  center-of-pressure  position  tends  toward  the 
centroid  especially  at  very  high  angles  of  attack.  If  a  value 
of  x/c  for  a  wing  at  some  large  angle  of  attack  such  as  45° 
were  known,  then  the  values  given  by  the  above  method  for 
the  range  of  shock  attachment  could  be  faired  to  the  a  =  45° 
value  to  obtain  an  estimate  over  a  broad  angle-of-attack  range. 
With  this  approach  in  mind,  the  center-of-pressure  positions 
for  a  =  45°  have  been  plotted  versus  Mach  number  for  £  1.6 
in  Figure  34  for  the  wings  of  the  Fidler-Baker  data  base 
(ref.  5)  and  the  Stallings-Lamb  data  base  (ref.  3).  Curves 
have  been  faired  through  the  data  of  the  Stallings-Lamb  data 
base  for  design  purposes.  The  Stallings-Lamb  wings  were 
chosen  because  the  Fidler-Baker  wings  were  tested  as  semispan 
models  on  a  reflection  plane  and  may  be  affected  by  boundary- 
layer  separation  at  high  angles  of  attack.  In  any  event, 
there  is  not  much  difference  between  the  two  sets  of  data  in 
Figure  34.  It  is  noted  that  the  effect  of  Mro  is  small  in  the 
range  of  1.6  to  4.6  shown  in  the  figure.  For  the  X  =  1 
wings  the  center  of  pressure  is  somewhat  forward  of  x/cr  =  0.5, 
the  centroid  position.  For  the  X  =  0  wings  it  is  also  forward 
of  0.667,  the  centroidal  position. 

Consider  now  results  for  J?/cr  at  a  =  45°  for  Mw  £  1.6 
as  shown  in  Figure  35.  The  results  for  flR  =  0.5  shown  in 
Figure  35(a)  show  definite  effects  of  Mach  number  on  x/cr 
below  a  Mach  number  of  unity.  An  effect  of  wing  section  on 
x/cr  at  a  =  45°  in  the  lower  transonic  range,  M^  =  0.6  to  1.0, 
is  evident  for  the  X  =  0.5  wing.  For  7R  = 1  the  wings  of 
Figure  35(b)  show  a  significant  effect  of  wing  section  for  the 
X  =  1  wings  only  (around  M^  -  1.0).  The  same  is  true  for  the 
/R  = 2  wings  data  shown  in  Figure  35(c).  Whether  this  effect 
is  a  true  transonic  effect  or  is  due  to  problems  of  testing 
a  semispan  model  on  a  reflection  plate  is  not  known. 


For  help  in  estimating  the  x/cr  positions  for  a  wing 
which  has  an  angle-of-attack  region  of  attached  flow  at  the 
leading  edge  in  the  supersonic  range,  the  data  in  Figures  34 
and  35  is  presented  in  tabular  form  in  Table  1.  This  table 
includes  the  3*  =  4  results  not  given  in  the  figures  but  which 
are  needed  for  wings  with  3*  >  2 .  The  data  are  sufficient  to 
cover  the  range  0.5  £  3*  <4.0  and  0  <  A  <  1.0.  For  wings  which 
do  not  fall  on  the  precise  values  of  31  and  A  of  the  wings  in 
the  data  base,  interpolation  is  to  be  used.  Generally, 
interpolation  is  best  done  linearly  in  taper  ratio,  inversely 
as  aspect  ratio,  and  inversely  as  /m^  -  1.  Normal-force  curves 
for  the  Fidler-Baker  data  base  are  available  in  reference  5, 
and  normal-force  curves  for  the  Stalling-Lamb  wings  are 
available  for  M  >  1.6  in  reference  3.  Tabulated  normal- 

CO  - 

force  coefficients  for  the  10  Stalling-Lamb  wings  are  given 
in  Table  2  for  M  >1.6. 

CO  — 

4.2.3  Detached  shock  case.-  The  detached  shock  case 
covers  all  wings  for  which  the  shock  is  not  attached  to  the 
leading  edge  for  any  angle  of  attack  when  considering  two- 
dimensional  flow  in  streamwise  planes  normal  to  the  wing  plan- 
form  and  parallel  to  the  root  chord.  The  leading  edge 
semiapex  angle  for  which  the  shock  wave  detaches  at  a  =  0°  is 
given  in  Figure  36. 

We  have  seen  that  wing  thickness  effect  can  have  a  large 
influence  on  the  axial  center-of-pressure  position  when  the 
leading-edge  shock  is  attached.  When  the  leading-edge  shock 
is  not  attached  even  at  a  =  0°  the  data  at  Moo  =  1.6  in  Figures  26 
and  29-33  still  exhibit  large  effects  of  wing  thickness  (airfoil 
section)  on  axial  center-of-pressure  position.  Let  us  now 
inquire  how  we  might  determine  this  quantity  theoretically. 

Linear  theories  based  on  the  wave  equation  for  determin¬ 
ing  the  center-of-pressure  position  of  a  wing  alone  with  or 
without  thickness  at  supersonic  speed  are  well  known  . 


However,  we  know  that  the  use  of  linear  theory  to  obtain  the 

velocity  field  combined  with  Bernoulli's  equation  to  calculate  j 

the  pressure  does  not  accurately  account  for  thickness  effects 

on  (x/cr)w  for  attached  leading-edge  shocks.  There  is 

thus  no  promise  this  procedure  will  do  so  for  detached  shocks. 

In  fact,  no  proven  method  is  now  available  for  determining 
thickness  effects  on  x/cr  for  detached  shocks.  A  method 
based  on  the  Euler  equations  has  promise  but  is  not  currently 
available.  Thus,  we  must  rely  for  the  present  on  experiment. 

With  respect  to  the  experimental  approach,  there  seems 
only  a  few  alternatives  available.  First  one  might  use  a  fin  ; 

for  which  the  wing  alone  corresponds  closely  to  one  of  the 
wings  of  the  Fidler-Baker  or  Stallings-Lamb  data  bases  or  may 
be  found  by  interpolating  within  them.  Again  the  fin  may  * 

correspond  to  a  wing  alone  for  which  data  are  otherwise 

available.  As  a  final  resort,  one  could  measure  the  properties  - 

of  the  wing  alone  in  a  wind  tunnel. 

8 

4.3  Effect  of  Vortices  on  Fin  Axial  Center- 
of-Pressure  Position 

* 

4.3.1  Panel  method  approach.-  Body  or  canard  vortices 
can  have  a  significant  effect  on  the  hinge  moment  of  an  ? 

all-movable  control  if  they  pass  in  close  proximity  to  it. 

This  can  occur  when  a  fin  is  rolled  well  up  into  the  first  or  «; 

second  quadrants  at  the  moderate  to  large  angles  of  attack  for 
which  body  vortices  exist.  It  can  also  occur  when  the  vortices  - 

from  canard  fins  pass  near  an  empennage.  It  is  not  known  how  ' 

to  determine  the  effect  of  such  vortices  on  axial  center-of- 
pressure  position  using  reverse -flow  methods  as  is  possible  for  / 

normal  force  and  lateral  center-of-pressure  position.  We  are 
therefore  interested  in  developing  a  simple  means  for  determin¬ 
ing  the  effect  of  vortices  on  axial  center-of-pressure 
positions. 

(• 
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In  order  to  gain  the  necessary  insight  to  model  the 
phenomenon  in  question,  we  will  make  calculations  using  a  panel 
method  (ref.  12)  to  study  the  vortex  effects.  For  this  purpose 
we  have  chosen  the  configuration  shown  in  Figure  37  for 
which  fin  data  are  available  in  reference  37.  In  the  figure 
comparison  is  made  between  data  and  the  prediction  of  the 
panel  method  for  fin  normal-force  coefficient,  root  bending- 
moment  coefficient,  and  hinge-moment  coefficient.  The 
comparisons,  for  =  1.6  for  fin  deflections  of  0°  and  -15°, 
show  excellent  agreement  for  fin  normal  force  and  bending 
moment  but  not  for  hinge  moment.  The  4>f=45°  results 
include  significant  body  vortex  effects,  and  the  <J>^  =135° 
results  are  free  of  such  effects.  Nevertheless  the  hinge- 
moment  results  for  the  <t>f  =135°  fins  are  not  well  predicted. 
This  result  is  associated  with  the  fact  that  panel  methods  do 
not  properly  account  for  the  second-order  effects  of  thickness 
on  axial  center-of-pressure  locations  as  discussed  previously. 
(This  configuration  will  be  discussed  subsequently  with  such 
second-order  effects  accounted  for.)  We  feel,  however,  that 
the  increment  in  wing  loading  due  to  the  body  vortices  is 
probably  well  predicted  by  the  panel  method. 

In  order  to  obtain  accurate  axial  center-of-pressure 
positions  one  might  inquire  how  many  chordwise  panels  are 
required.  A  systematic  set  of  calculations  has  been  made 
varying  the  number  of  chordwise  panels  from  3  to  7,  while 
keeping  the  number  of  spanwise  rows  equal  to  5.  The  calcula¬ 
tions  were  made  for  the  example  configuration  of  the  previous 
figure  for  the  leeward  fin  at  =45°.  Calculations  were  made 
with  and  without  body  vortices.  The  vortex  strengths  and 
positions  were  estimated  from  correlation  curves  in  refer¬ 
ence  17.  Figure  38  presents  the  results  of  the  calculations 
for  the  center-of-pressure  positions  with  and  without  body 
vortices  plotted  versus  1/n,  where  n  is  the  number  of  chordwise 
panels.  It  is  noted  that  without  vortices  the  value  of  x/cr 


is  close  to  that  obtained  from  DATCOM  for  all  values  of  1/n. 

The  body  vortices  cause  about  a  3.5  percent  forward  shift  in 
the  fin  center-of-pressure  position  for  all  values  of  1/n, 
a  shift  which  can  signficantly  influence  hinge  moment. 

A  series  of  calculations  has  been  made  for  the  same  con¬ 
ditions  as  those  of  Figure  38,  but  the  fin  taper  ratio  has  been 
varied  from  A  =  0  to  0.423  to  1.0  with  a  fixed  exposed  semi¬ 
span  of  11.43  as  shown  in  Figure  37.  The  results  plotted  in 
Figure  39  form  the  basis  for  the  simplified  method  of  deter¬ 
mining  vortex  effects  on  fin  axial  center-of-pressure  position. 
Figure  39(a)  for  A  =  0  shows  the  local  center-of-pressure 
positions  at  five  spanwise  stations  with  and  without  body 
vortex  effects.  The  lines  joining  these  positions  are  nearly 
straight  and  can  be  thought  of  as  lifting  lines.  Consider  now 
the  wing  center  of  pressure  without  vortices  and  the  center  of 
pressure  of  the  vortex  normal  force.  The  vortex  normal-force 
center  of  pressure  is  inboard  and  lies  close  to  the  lifting 
line  without  vortex  effects.  Since  the  vortex  induced  normal 
force  is  downward,  the  fin  center  of  pressure  is  moved  aft 
by  vortex  effects. 

Examine  now  figure  39(b)  for  a  fin  with  A  =  0.423.  The 
results  are  qualitatively  similar  to  those  for  A  =  0,  and  a 
significant  shift  of  the  fin  center  of  pressure  in  a  rearward 
direction  is  due  to  the  body  vortices.  Examine  now  Figure  39(c) 
for  the  rectangular  fin  with  the  same  exposed  semispan  as  the 
other  fins.  The  lifting  lines  are  now  essentially  unswept, 
and  the  lifting  line  for  no  body  vortex  effects  is  not  straight. 
The  center  of  pressure  of  the  vortex  normal  force  lies  slightly 
behind  this  lifting  line,  unlike  the  situation  for  the  other 
two  cases.  However,  because  the  lifting  lines  are  nearly 
unswept,  the  effect  of  the  vortices  on  axial  center-of-pressure 
location  is  only  about  1  percent  of  the  root  chord. 


4.3,2  Approximate  engineering  method.-  The  foregoing 
results  suggest  assumptions  that  lead  to  a  simple  method 
for  determining  the  effect  of  the  vortices  on  fin  center-of- 
pressure  positions.  It  is  first  assumed  that  the  lifting 
line  is  a  straight  line  passing  through  the  fin  center-of- 
pressure  position  without  vortex  effects  at  a  constant 
percent  local  chord.  It  is  also  assumed  that  the  center  of 
pressure  of  the  vortex  normal  force  lies  on  this  lifting  line, 
We  can  get  an  idea  of  the  degree  of  approximation  of  this 
assumption  from  the  following  values  for  three  cases  of 
Figure  39  for  the  center-of-pressure  locations  in  terms  of 
percent  local  chord. 


0.423 

1.00 


(x-x£e)/c 

No  Vortices 

0.452 

0.459 

0.439 


(xv“x£e)/c 


0.423 


0.427 

0.491 


c  =  local  chord  through  center  of  pressure 

We  have  values  for  the  coordinates  of  x,y  of  the  fin  center- 
of-pressure  position  with  no  vortex  effect  from  the  computa¬ 
tional  method.  We  obtain  yv  due  to  the  loading  induced  on 
the  fin  by  the  vortex  (vortices)  in  the  program  using  reverse 
flow  methods.  Knowing  the  value  of  (xv  -  x£e^c  the  ai3°ve 
assumption  and  knowing  yv,  we  have  fixed  xv. 

We  now  list  the  steps  required  to  determine  the  center  of 
pressure  of  the  vortex  induced  normal  force  as  a  fraction 
of  the  root  chord  (x/cr)  . 


center-o 
pressure 
vortex  C 


1.  Determine  the  lateral  center-of-pressure  position 
of  the  vortex  normal  force,  (y/a)v,  from  MISSILE  2  as  now 
constituted. 


2.  Determine  (x/cr)w^  by  the  methods  discussed  in 
the  previous  section  (section  4.2)  for  C^p  with  no  vortices. 


3.  Determine  the  lateral  center  of  pressure  of  the  fin 
(with  or  without  deflection)  and  with  no  vortex  effects 

(y  -  a)/(sm  -  a)w^  by  interpolating  or  extrapolating  in  the 
results  of  Program  MISSILE. 

4.  Determine  the  local  fractional  chord  t  of  the  fin 
center  of  pressure  with  no  vortex  effects  from 


5.  Determine  the  percent  root  chord  corresponding  to  the 
axial  position  of  the  normal  force  induced  on  the  fin  by  the 
vortices,  (x/cr)v. 


=  t  +  (1  -  t)  (1  -  A) 


sm  -  a 


We  know  the  normal  force  induced  on  the  fin  by  the 
vortices  and  we  are  now  in  a  position  to  determine  its 
contribution  to  the  hinge-moment  coefficient.  It  is 
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noted  that  in  this  approximate  method  the  normal  force  and 
hinge  moment  induced  by  the  body  vortices  go  to  zero  together. 

No  induced  camber  effects  can  be  handled. 

4.4  Other  Changes  in  Program  MISSILE 

4.4.1  Introductory  remarks. -  We  have  made  a  number  of 
changes  in  MISSILE  2  for  more  accurate  calculation  of  hinge 
moments.  The  first  major  change  is  that  we  now  require  as 
input  a  table  which  gives  the  wing-alone  axial  center-of- 
pressure  position  as  a  function  of  fin  normal-force  coefficient 
as  well  as  the  normal-force  coefficient  as  a  function  of 
angle  of  attack.  Methods  of  interpolating  or  extrapolating 

in  the  wing-alone  data  base  to  obtain  this  information  for  a 
wing  not  in  the  data  base  are  given  below.  Methods  for 
extrapolating  the  fin  normal  force  with  no  vortex  effects  as 
given  in  the  data  base  to  high  aspect  ratios  and  high  Mach 
numbers  have  been  improved.  Also,  the  application  of  the 
equivalent  angle  of  attack  concept*  has  been  modified  to  account 
for  vortex  effects  prior  to  accounting  for  fin  deflection 
effects.  The  two  effects  are  commutative  only  for  a  linear 
normal-force  curve.  Finally,  improvements  have  been  made  in 
the  treatment  of  fin-fin  interference  due  to  control  deflection 
using  the  equivalent  angle  of  attack  concept  based  on  certain 
physical  arguments. 

4.4.2  Wing-alone  characteristics  as  input  data. - 

1.  Effect  of  Mach  number  on  normal  force.  In  case  we 
do  not  have  data  for  the  wing  alone  at  the  conditions  of 
interest,  we  may  choose  to  interpolate  (or  extrapolate)  in 
the  Fidler-Baker  or  the  Stallings-Lamb  data  bases.  Such 
interpolations  may  be  made  with  respect  to  JR,  X,  or  M^.  With 
regard  to  extrapolating  above  Mra  =  3.0  in  the  Fidler-Baker  data 

_ 

See  Appendix  C  for  a  description  of  the  equivalent  angle- 
of-attack  concept. 


base,  the  Stallings-Lamb  data  base  which  goes  up  to  Mffi  =  4.6 
sheds  light  on  this  question.  In  Figure  40  the  wing  normal- 
force  coefficient  from  Table  2  has  been  plotted  versus  Mm  for 
constant  angle  of  attack  for  wings  with  ®=0.5,  A  =  l,  and 
i“R  =4.0,  A  =  0.5  as  extreme  cases.  Straight  line  projections 
have  been  drawn  through  the  =  2.16  and  2.86  data  points 
to  illustrate  the  accuracy  of  straight-line  extrapolation  to 
higher  Mach  numbers.  It  is  clear  that  significant  errors  can 
be  made  if  such  extrapolation  is  carried  out  very  far  above 
M  =  3.0. 

oo 

A  better  extrapolation  means  should  take  advantage  of  the 

known  fact  that  the  aerodynamic  coefficients  change  slowly 

with  Mach  number  and  become  asymptotic  as  M^  -*  °°.  A  parameter 

2  -1/2  °° 

which  has  this  property  is  (M  -  1)  .  The  data  of  Figure  40 

2  °°  -1/2 

have  been  plotted  versus  (M^-l)  '  in  Figure  41.  Extrapolat¬ 

ing  to  higher  Mach  numbers  based  on  the  M^  -•  2.16  and  2.86  data 
shows  great  improvement  in  accuracy  over  the  linear  extrapola¬ 
tion.  As  a  consequence  we  have  used  this  type  of  interpolation 
in  the  computer  program  for  data  above  M^  =  3  for  both  the 
wing  alone  and  the  fin  in  the  presence  of  the  body.  The 
linear  interpolation  is  retained  within  the  Mach  number  range 
of  the  fin-on-body  data  base,  0.8  £  M^  £  3.0. 


If  the  Mach  number,  aspect  ratio,  and  taper  ratio  do  not 
correspond  to  values  in  the  wing-alone  data  base,  it  is 
suggested  that  linear  interpolation  be  performed  in  the  data 
base.  The  three  closest  wings  to  the  subject  wing  are  selected 
and  their  characteristics  are  interpolated  to  the  correct  Mach 
number  before  the  aspect-ratio,  taper-ratio  interpolation. 

If  the  Mach  number  of  interest  is  outside  the  realm  of  the 
data  base,  that  Mach  number  extrapolation  is  done  as  described 
above  before  aspect-ratio,  taper-ratio  interpolation. 

2 .  Effect  of  Mach  number  on  axial  center-of-pressure 
location.  A  study  similar  to  that  for  extrapolation  of  normal 
force  has  been  carried  out  for  extrapolation  of  axial 


center-of-pressure  location.  The  results  are  shown  in  Figures 
42  and  43.  They  show  that  axial  center-of-pressure  positions 
cannot  be  extrapolated  in  Mach  numbers  above  3  with  confidence 
either  linearly  or  as  1//m£  -  1.  Generally,  the  methods  of 
predicting  x/cr  previously  described  should  be  relied  on. 

For  Mot  >  3  practically  all  wings  will  have  a  region  of  attached 
leading-edge  shock  waves. 

3 .  Extrapolation  of  normal-force  coefficient  for  /R  >2.0. 
There  is  a  requirement  to  extrapolate  the  wing  alone  and  fin- 
on-body  normal-force  coefficients  to  aspect  ratios  greater 
than  2.  Since,  for  a  given  angle  of  attack,  normal-force 
coefficient  becomes  asymptotic  to  the  two-dimensional  value 
as  /R  -*■<»,  it  is  not  accurate  to  extrapolate  linearly  in  /R 
without  limit.  Instead,  extrapolation  in  1/®  has  the  correct 
asymptotic  behavior.  Consider  that  data  are  available  for 
iR  =  1  and  /R  =  2  and  we  desire  the  extrapolation  to  iPR  >  2 .  We 
can  use  the  Stallings-Lamb  data  of  Table  2  to  test  our 
extrapolation  scheme  for  X  =  0.5  since  data  are  available  to 
®  =4.  Wing-alone  normal-force  coefficient  is  plotted  versus 
1/iR  for  constant  a  in  Figure  44  for  two  Mach  numbers  for 
X  =  0.5  wings.  It  is  seen  that  straight  line  extrapolation 
predicts  the  iR  =  4  results  quite  well  in  all  cases. 


for  jR  >  2.  Plots  were  made  similar  to  Figure  44  to  test  the 
feasibility  of  extrapolating  axial  center-of-pressure  results 
to  aspect  ratios  greater  than  iR  >  2 .  The  Stallings-Lamb  wing 
data  for  X  =  0.5  and  M  =1.6  and  M  =3.5  are  plotted  versus 
1//R  in  Figures  45(a)  and  45(b).  The  lines  connecting  the 
/R  =  1  and  JR  =  2  data  extended  to  &  =  4  sometimes  provide  a 
good  prediction  and  sometimes  do  not.  Accordingly,  we  do 
not  recommend  extrapolating  axial  center  of  pressure  in  aspect 
ratio. 
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5.  Interpolation  of  wing-alone  Cn  and  x/cr  in  wing 
taper  ratio.  The  question  of  interpolating  in  a  data  base  to 
obtain  wing  alone  or  fin-on-body  CN  and  x/cr  for  values  of 
X  not  in  the  data  base  arises  in  preparing  input  data  for  the 
program  and  in  the  program  itself  in  obtaining  fin-on-body 
normal  force.  Generally,  values  of  X  run  from  0  to  1.0,  and 
the  question  of  extrapolating  to  values  greater  than  1  does 
not  arise.  Data  exist  for  the  Stallings-Lamb  wings  to  show 
the  effect  of  X  on  the  measured  aerodynamic  quantities  for 
PR  =0.5,  1.0,  and  2.0.  We  show  the  effect  of  X  on  wing  CN 
for  ®  =  2  wings  in  Figure  46  for  =  1.6  and  Mot  =  3.50.  The 
curves  drawn  through  the  data  have  very  little  curvature  so 
that  straight  line  interpolation  between  data  points  is 
sufficiently  accurate  for  our  purposes. 

Considering  now  axial  center-of-pressure  position, 
this  quantity  is  shown  for  the  same  parameters  in 
Figures  47(a)  and  47(b).  Curves  are  faired  through  the  data. 
Straight  lines  between  the  data  points  deviate  from  the  faired 
curves  by  less  than  0.01.  Thus,  linear  interpolation  in  X  for 
x/cr  should  be  adequate  for  hinge-moment  purposes  in  most 
cases. 


4.4.3  Interpolation  and  extrapolation  for  the  fin-on- 
body  data  base.-  The  fin-on-body  data  base  is  needed  solely 
for  finding  the  normal  force  of  the  fin  on  the  body  including 
all  interference  effects  except  body  vortices;  it  is  not  used 
to  determine  the  axial  center-of-pressure  location  of  the  fin 
on  the  body.  Accordingly,  we  are  interested  only  in  inter¬ 
polation  or  extrapolation  with  respect  to  fin  normal-force 
coefficient.  Much  of  section  4.4.2  applies  to  this  question 
because  under  the  equivalent  angle-of-attack  concept  a 
parallelism  has  been  set  up  between  the  nonlinear  character¬ 
istics  of  the  wing  alone  and  the  fin  in  the  presence  of  the 
body. 
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Consider  the  aspect  ratio,  Mach  number  domain  of  the 
present  data  base  shown  in  Figure  48.  The  aspect  ratio- 
Mach  number  combinations  for  which  data  are  available  in  the 
present  data  base  are  shown  as  circles,  and  the  domain  of  the 
present  data  base  is  contained  within  the  boundary  enclosing 
these  circles.  The  dotted  boundary  is  that  for  a  tri-service 
effort  presently  underway  which  can  be  used  to  extend  the 
domain  of  application  of  the  present  effort  at  a  later 
date  (ref.  38) . 


With  regard  to  interpolation  within  the  present  data 
base,  we  retain  linear  interpolation  in  Mra,  ®  ,  and  X  as 
in  MISSILE  2.  However,  the  new  feature  added  to  the  hinge- 
moment  program  is  the  ability  to  extrapolate  in  /R  and  M^  in 
the  directions  shown  in  the  figure  by  a  more  accurate  means. 
The  extrapolation  in  aspect  ratio  is  done  inversely  on  A* . 


The  extrapolation  with  respect  to  M^  has  been  done  proportional 
to  (M2  -  1)"1/2. 


Another  change  has  been  introduced  into  the  hinge-moment 
program.  In  Program  MISSILE  an  approximate  wing-alone  CN 
data  base  has  been  built  into  the  program.  This  together  with 
the  fin-on-body  CN  data  base  gives  values  of  Kw  used  in  the 
equivalent  angle-of-attack  method.  Since  for  the  hinge-moment 
program  the  wing-alone  normal-force  curve  is  input  data,  it  is 
necessary  to  reevaluate  Kw.  This  is  done  internal  to  the 
program.  Since  wing-alone  normal  force  is  not  influenced 
greatly  by  airfoil  section,  this  change  should  not  influence 
the  calculated  hinge  moment  very  much  (except  at  transonic 
stall)  and  is  done  for  internal  consistency. 


The  present  model  for  the  strengths  and  positions  of  the 
body  vortices  includes  Mach  number  up  to  3.0.  For  higher 
Mach  numbers  linear  extrapolation  is  used  based  on  the  M^  =  2 

and  M  = 

00 


3  results. 


4.4.4  Changes  in  implementation  of  equivalent  angle- 


of-attack  concept.-  We  will  now  describe  two  changes  in  the 
implementation  of  the  equivalent  angle-of-attack  concept  which 
have  been  incorporated  into  the  hinge-moment  program  and 
which  are  not  in  MISSILE  2.  A  detailed  description  of  the 
equivalent  angle  of  attack  with  instructions  for  its  use 
are  given  in  reference  38. 

Briefly  the  equivalent  angle  of  attack  is  the  angle  of 
attack  of  the  wing  alone  for  which  it  has  twice  the  normal 
force  of  the  fin  in  the  presence  of  the  body.  The  equivalent 
angle  of  attack  of  the  fin  depends  on  body  upwash  (angle  of 
attack  and  roll  angle) ,  body  vortex  induced  downwash  on  the 
fin,  and  wash  induced  on  the  fin  by  deflection  of  other  fins. 
The  equivalent  angle  of  attack  depends  on  the  order  of 
application  of  the  vortex  effects  and  fin  deflection  effects 
because  the  fin-fin  interference  factors  entering  the  method 
depend  on  angle  of  attack.  In  Program  MISSILE  the  body 
vortex  effects  were  imposed  after  the  fin  deflection  effects. 

In  this  hinge-moment  program  this  procedure  has  been  reversed. 
The  fin-fin  interference  factors  for  use  in  the  method  can 
then  be  determined  from  wind-tunnel  data  without  first 
removing  body  vortex  effects  from  the  data,  a  process  of  some 
uncertainty. 

The  second  change  made  in  anticipation  of  the  tri-service 
data  base  is  simplification  of  the  fin-fin  interference  treat¬ 
ment.  In  Program  MISSILE  the  effect  of  panel  deflections  on 
the  equivalent  angle  of  attack  of  a  particular  panel  was 
obtained  by  a  tangent  addition  theorem  which  did  not  permit 
the  equivalent  angle  of  attack  to  exceed  90°.  However,  it 
is  clear  that  fin  deflection  increases  the  fin  equivalent 
angle  of  attack  linearly  and  that  an  equivalent  angle  of  attack 
greater  than  90°  is  possible.  Thus,  if  aecj^  equivalent 

angle  of  attack  of  fin  i  with  no  fin  deflection  and  aeq^  is 
that  with  all  fins  deflected,  then 
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where  A^j  are  fin-fin  interference  factors.  For  instance, 

A^j  is  the  amount  the  equivalent  angle  of  attack  of  fin  i  is 
increased  by  unit  deflection  of  fin  j.  The  factor  A^  should 
be  close  to  unity.  Theoretical  values  of  A^^  obtained  from 
slender-body  theory  (ref.  1)  are  now  used  in  the  hinge- 
moment  program  with  corrections  for  how  much  of  fin  i  lies 
within  the  zone  of  influence  of  fin  j.  The  zone  of  influence 
correction  is  also  used  in  the  present  version  of 
MISSILE  2. 


5.  FINS  OF  GENERAL  PLANFORM 

As  a  result  of  the  present  study,  a  method  can  be 
specified  for  application  to  fins  of  general  planform.  The 
method  will  be  described  for  possible  future  developement  and 
verification. 

By  a  fin  of  general  planform  we  mean  one  for  which  the 
leading,  trailing,  and  side  edges  can  be  approximated  by 
straight  line  segments.  The  airfoil  sections  should  be  such 
that  their  contours  can  be  approximated  by  straight  line 
segments.  Its  leading  edges  should  be  sharp. 

The  method  is  meant  to  apply  only  to  the  range  of  leading- 
edge  shock  attachment  in  a  -  Mro  space.  The  criterion  for  shock 
attachment  shall  be  applicability  of  shock-expansion  theory 
in  planes  normal  to  the  fins  and  parallel  to  the  body  axis. 

The  first  task  in  the  method  is  to  determine  the  normal- 
force  coefficient  and  axial  center  of  pressure  for  the  wing 
alone  using  the  panel  theory  of  reference  12  with  no  thickness 
effects  in  the  calculation.  These  quantities  are  determined 
up  to  an  angle  of  attack  of  about  20°.  Some  nonlinearities 
are  accounted  for  by  the  method.  Then  the  shock-expansion 
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theory  computer  program  to  determine  the  shift  in  axial 
center-of-pressure  position  due  to  wing  thickness  is  applied. 
This  establishes  the  axial  center-of-pressure  position  up  to 
leading-edge  shock  detachment.  An  estimate  of  x/cr  at  a  =  45° 
is  made  next.  This  is  done  by  determining  the  wing  planform 
centroid,  and  then  subtracting  an  estimated  quantity  from  it 
using  the  results  of  Figures  34  and  35  as  a  guide.  A  hand 
fairing  between  the  results  of  shock-expansion  theory  and  that 
for  a  =  45°  is  then  made. 

The  second  main  task  is  to  choose  an  equivalent  wing  in 
the  fin-on-body  data  base  to  replace  the  general  planform  in 
determining  the  normal  force  on  the  fin  in  the  presence  of  the 
body.  (If  the  angle  of  attack  range  is  below  20°,  we  can 
circumvent  this  step  by  using  the  panel  program  for  the  wing- 
body  combination  to  calculate  this  quantity  directly.)  However, 
if  we  desire  to  use  the  engineering  method,  we  must  choose  an 
equivalent  wing  for  normal  force  in  the  following  way. 

1.  The  normal  force  of  the  equivalent  clipped  delta 

wing  should  equal  that  of  the  wing  of  general  plan- 
form  at  the  same  angle  of  attack  and  Mach  number. 


S  C„  —  S,CM 
e  N  f  N- 

e  t 

S  ~  reference  area 

e  ~  equivalent  wing 

f  ~  wing  of  general  planform 

2.  The  exposed  semispans  shall  be  equal. 

The  first  requirement  is  based  on  the  assumption  that  if 
the  wing-alone  normal  forces  have  the  same  magnitude,  then  the 
normal  forces  acting  on  the  corresponding  fins  on  the  body  will 
be  the  same.  The  assumption  is  suggested  by  the  following 
arguments.  The  normal  force  on  a  fin  in  the  presence  of  the 
body  is  related  to  that  on  the  fin  as  part  of  the  wing  alone 


by  a  factor  Kw  which  depends  principally  on  the  a/sm  ratio 
independent  of  planform.  This  condition  is  met  by  having 
the  fin  of  general  planform  and  the  equivalent  fin  have  the 
same  exposed  semispans.  Also  the  effect  of  body  vortices  on 
the  fin  normal  force  depends  on  a  reverse  flow  method  which 
uses  weighting  factors  depending  only  on  a/sm  and  which  are 
independent  of  planform.  In  applying  the  first  rule  it  may 
be  enough  to  match  normal-force  curve  slope  times  reference 
area.  The  whole  purpose  of  the  "equivalent  wing"  is  to  use 
the  fin-on-body  data  base  in  the  program  to  get  an  estimate  of 
the  normal  force  of  the  fin  of  general  planform  mounted  on  the 
body.  Any  considerations  of  planform  on  axial  center-of- 
pressure  position  are  extraneous  since  this  is  handled  through 
the  wing-alone  data  which  are  input. 

The  third  main  task  in  the  procedure  is  to  apply  MISSILE  2 
as  modified  to  the  problem  using  the  equivalent  wing  to  get  fin 
normal-force  coefficient.  This  quantity  and  the  axial  center-of 
pressure  position  estimated  for  the  fin  of  general  planform  then 
yield  the  hinge-moment  coefficient.  In  case  the  fin-on-body 
normal  forces  are  calculated  by  the  panel  program,  the  axial 
center-of-pressure  position  estimated  for  the  fin  of  general 
planform  should  still  be  used  in  calculating  the  hinge  moment 
rather  than  that  from  the  panel  program.  In  this  step  a 
refinement  may  be  necessary.  It  may  be  necessary  to  determine 
the  incremental  distance  between  the  axial  center-of-pressure 
positions  of  the  vortex  normal  force  and  normal  force  without 
vortices  present  by  using  the  panel  method.  The  normal  force 
for  no  vortices  would  then  be  applied  at  the  axial  center-of- 
pressure  position  estimated  for  the  fin  of  general  planform 
and  the  vortex  normal  force  would  be  applied  at  the  incremental 
distance  from  this  position  determined  on  the  basis  of  the 
panel  method. 


At  the  present  time  it  is  not  clear  how  to  treat  the  fin 
of  general  planform  in  the  case  where  the  shock  is  detached 
from  the  leading  edge  at  zero  angle  of  attack  unless  exper- 
mental  wing-alone  data  on  normal  force  and  axial  center-of- 
pressure  positions  are  available. 

6.  COMPARISON  OF  EXPERIMENTS  AND  THEORY 
6.1  Introductory  Remarks 

Comparisons  will  be  made  between  prediction  and  experi¬ 
ment  for  a  number  of  cases  to  investigate  the  accuracy  and 
range  of  applicability  of  the  method.  The  range  of  applica¬ 
bility  will  be  defined  in  terms  of  wing  aspect  ratio,  wing  taper 
ratio,  angle  of  attack,  fin  deflection,  and  Mach  number.  In 
addition,  roll  angle  is  a  parameter.  The  configurations  will 
be  principally  cruciform  wing-body  combinations  although  we 
also  consider  an  exploratory  application  to  a  wing-body- tail 
combination. 

Two  sets  of  data  are  particularly  useful  for  making  the 
comparisons.  The  first  set  is  that  of  Lamb  and  Trescot 
(ref.  39)  and  of  Lamb,  Sawyer,  Wassum,  and  Babb  (refs.  37, 

40,  and  41).  These  data  are  particularly  interesting  because 
the  fin  pressure  distributions  are  available  to  diagnose  any 
significant  differences  between  experiment  and  theory.  The 
configurations  investigated  are  shown  in  Figure  49.  The  second 
set  of  data  is  that  for  the  Army  Generalized  Missiles  described 
in  reference  42  and  shown  in  Figures  50  and  51.  While  no 
pressure  distribution  data  are  available,  there  are  fin  force 
and  moment  data  for  wing-body  and  wing-body-tail  combinations 
for  several  different  fins  at  a  number  of  roll  angles  and  at 
subsonic  and  supersonic  speeds.  Other  data  exist  in  the  public 
domain,  but  they  are  not  so  extensive  nor  so  systematic  as  the 
foregoing  sources.  We  start  by  considering  the  first  set  of 
data. 


6.2  Data  of  Lamb  and  Trescot 

6.2.1  Wing-alone  characteristics.-  Comparisons  will  be 
shown  for  fin  B  of  Figure  49  at  Mach  numbers  of  1.6,  2.35, 
and  3.7.  The  first  task  is  to  find  the  wing-alone  values 
of  CN  and  x/cr  for  this  fin  at  the  various  Mach  numbers  for 
input  into  the  program.  To  obtain  the  normal-force  character¬ 
istics,  we  have  used  interpolation  in  the  Stallings-Lamb 
data  base  (ref.  3)  .  The  geometric  parameters  for  fin  Tg 
are: 

m  ~  2.812 
X  =  ct/cr  =  0.423 

a/sm  =0.25 

xttt/cr  =  °*589 

The  normal-force  curves  of  the  Stallings-Lamb  data  base 
have  been  fitted  by  least  squares  with  a  truncated  sine  series 
of  the  following  form. 

CN  =  A-^sina  +  A2sin3a  +  A3sin5a 

The  values  of  the  coefficients  A^,  A2 ,  and  A^  are  given  in 
Table  3.  Good  fits  to  the  data  were  obtained  in  all  cases. 

To  obtain  the  values  of  A^ ,  A2 ,  A^  for  TB,  linear  interpola¬ 
tion  in  aspect  ratio  between  the  7R  =2,  X  =  0.5,  and  2R  =  4 
and  X  =  0.5  wings  was  used.  The  taper  ratio  effect  was 
obtained  by  linear  interpolation  between  the  7R  =2,  X  =  0  wing 
and  the  ®  =2,  X  =  0.5  wing.  Linear  interpolation  between 
Mach  numbers  was  used  for  the  particular  Mach  numbers  1.60, 
2.36,  and  3.70.  The  normal-force  curves  so  obtained  are  shown 
in  Figure  52. 

To  obtain  the  wing-alone  center-of-pressure  positions  for 
Tg,  interpolation  was  also  made  in  the  Stallings-Lamb  data  base 
As  has  been  previously  pointed  out,  the  wing-alone  axial  center 
of-pressure  position  is  influenced  by  airfoil  section. 
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Therefore,  the  dimensions  of  equivalent  wing  of  the  Stallings- 
Lamb  data  base  for  =2.812  and  X  =  0.423  were  found  by 
linear  interpolation.  The  semiapex  angle  normal  to  the  leading 
and  trailing  edges  of  the  data  base  wing  is  15°  while  for 
fin  Tg,  it  is  13.1°.  Since  the  wings  are  fairly  close,  it 
was  decided  not  to  correct  the  interpolation  data  for  center- 
of  pressure  positions.  Also,  shock-expansion  theory  will  not 
go  through  at  Mm  =  1.60  because  the  leading-edge  shock  is 
always  detached.  The  resulting  center-of -pressure  curves 
are  shown  in  Figure  53.  The  data  from  Figures  52  and  53  go 
into  the  computer  program  as  input. 

6.2.2  Comparison  of  experiment  and  theory.-  Comparisons 
have  been  made  between  prediction  and  data  for  fin  TB  for  the 
following  matrix  of  conditions: 

Mro  =  1.6,  2.36,  and  3.70 

0°  <  a  <  20° 

6  =  0° , -15° 

<pf  =  45°,  90°,  135° 

Consider  first  the  Mm  =  1.60  results  shown  in  Figure  54.  The 
results  for  the  horizontal  fin  position  (<j>f  =  90°)  as  shown 
in  Figure  54(a)  and  54(b)  show  good  agreement  for  both  fin 
normal-force  and  hinge-moment  coefficients  for  6=0  and  -15°. 
The  same  quantities  for  the  windward  fin  location  (<j>f  =  135°) 
are  shown  in  Figures  54(c)  and  54(d).  The  6=0°  results  are 
in  good  agreement,  but  the  6 = -15° results  are  not  good  in  the 
high  angle-of-attack  range.  This  effect  appears  to  be  the 
result  of  fin  unporting  in  the  wing-body  juncture.  We  will 
subsequently  discuss  this  phenomenon.  Consider  the  comparison 
for  the  fin  on  the  leeward  side  of  <p^  =  45°  shown  in 
Figures  54(e)  and  54(f).  For  6  =  0  the  normal-force  results  are 
nonlinear  and  in  fair  accord  with  experiment  suggesting  that 
the  vortex  effects  have  been  fairly  well  modeled  for  normal 
force.  For  6=0  the  hinge-moment  nonlinearities  are  somewhat 
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stronger  than  predicted.  For  6  =  -15°  both  CNF  and  CHM 
well  predicted  at  low  a  but  neither  is  well  predicted  at  high  a. 

We  now  consider  what  the  pressure  distributions  show 
concerning  the  behavior  of  fin  TB  for  <f>  =  135°  and  6  =  -15° 
as  shown  in  Figures  54(c)  and  54(d).  To  obtain  an  idea  of  the 
upwash  distribution  at  the  fin  location  with  no  fin  present, 
we  have  calculated  it  using  the  supersonic  linear  theory  for 
the  body  alone.  Body  vortex  effects  are  negligible  for 
<pf  =  135°.  The  resultant  local  angles  of  attack  with  respect 
to  the  fin  chord  plane  are  shown  in  Figure  55  for  6=0  and 
6  =-15°.  The  pressure  distributions  from  reference  37  are 
shown  in  Figure  56  for  the  upper  and  lower  surfaces  at  0.02 
and  0.22  of  the  exposed  semispan.  Inspecting  first  the  pres¬ 
sure  on  the  leading  edge  wedge  for  y/s  =  0.02,  we  note  that  the 
upper  surface  has  pressures  about  equal  to  those  of  the  lower 
surface  despite  a  calculated  upwash  angle  of  10°.  This  result 
corresponds  to  nearly  zero  local  angle  of  attack  rather  than 
the  predicted  angle.  Behind  the  x/c  =  0.2  location,  the  data 
show  a  more  normal  behavior  with  regard  to  the  magnitude  of 
the  pressure  distributions  on  the  upper  and  lower  surfaces. 
However,  they  do  not  show  uniform  pressure  on  the  flat 
surfaces.  Inspecting  the  pressure  distributions  for  the 
y/s  =  0.22  station  reveals  the  pressure  distributions  to  be 
similar  to  what  might  be  expected  for  the  airfoil  at  a  slight 
positive  angle  of  attack.  It  thus  appears  some  wing-body 
juncture  phenomenon  occurs  which  nearly  disappears  by 
y/s  =0.22.  It  is  suspected  this  is  a  local  gap  effect 
associated  with  unporting  and  possibly  separation  of  the 
boundary  layer.  However,  pressure  distributions  at  the  next 
higher  Mach  number,  Mm  =  2.36,  show  that  the  nonlinear  effect 
has  substantially  disappeared. 

Curves  directly  analogous  to  those  of  Figure  54  for 
M^,  =  1.6  are  plotted  in  Figure  57  for  =  2.36  and  in 
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Figure  58  for  =  3.7.  The  predictions  for  M^  = 
extrapolating  in  Mro  out  of  the  present  fin-on-body  data  base. 
We  now  look  at  the  various  cases  as  a  function  of  Mach  number 
to  see  how  the  nonlinearities  exhibited  at  M^  =  1.6  change 
with  Mach  number. 

Comparing  parts  (a)  of  Figures  54,  57,  and  58,  we  see 
that  the  prediction  of  CNF  for  all  Mach  numbers  are  fairly 
good  except  that  the  6  =  -15°  case  shows  more  download  in  the 
low  a  range  than  predicted.  The  difference  between  the  6=0 
and  6  =  -15°  normal  forces  are  associated  with  the  fin  inter¬ 
ference  factor  for  deflection,  kw,  which  is  known  to  vary  with 
a  and  6,  but  no  data  base  for  this  variation  is  yet  available 
It  is  planned  to  obtain  such  data  in  the  tri-service  study 
(ref.  38). 

Consider  now  parts  (b)  of  the  figures.  We  see  that 
hinge-moment  coefficients  are  fairly  well  predicted,  and 
that  the  hinge-moment  curves  are  more  nonlinear  than  the 
normal-force  curves  as  a  result  of  movements  in  x/cr  with 
changes  in  CNF« 

Consider  now  the  fin  rolled  to  the  windward  side  of  the 
body  with  <J)f  =  135°,  and  examine  the  normal-force  curves  of 
parts  (c)  of  the  figures.  The  fin  normal-force  coefficient 
at  M^  =  1.6  and  2.36  is  well  predicted  for  6=0  but  is 
underpredicted  for  M^  =  3.7.  This  probably  is  the  result  of 
the  large  extrapolation.  What  is  of  interest  is  that  the 
change  of  CNF  between  6=0  and  6  =  -15°  is  underpredicted  by 
the  calculation.  At  the  two  higher  Mach  numbers  15°  of 
deflection  produces  more  normal  force  than  15°  of  angle  of 
attack  with  6=0.  This  again  points  to  the  desirability  of 
building  a  kw  data  base  into  the  prediction  method  when  it 
becomes  available.  Another  method  of  improvement  based  on  the 
panel  method  will  be  subsequently  discussed. 


Consider  now  parts  (d)  of  the  figures  which  show  hinge- 
moment  comparisons  for  <p^  =  135°.  For  6  =  0  the  prediction 
of  hinge-moment  coefficient  is  within  expected  accuracy. 
However,  for  6  =  -15°  the  predictions  are  satisfactory  only 
at  low  angles  of  attack.  A  means  of  improving  the  accuracy 
of  prediction  for  the  higher  angle  of  attack  range  for  <p~  =  135° 
will  be  described  subsequently. 

Consider  now  parts  (e)  of  Figures  54,  57,  and  58  which 
show  the  predicted  and  measured  fin  normal-force  coefficients 
for  (J>£  =  45°.  The  =  1.6  comparisons  for  6  =  0  are  good, 
but  the  predictions  for  6  =  -15°  are  not  good  at  the  higher 
angles  of  attack.  The  =  2.36  and  3.7  predictions  are  not 
good  at  all.  The  reason  for  the  discrepancy  is  not  known. 

The  hinge-moment  results  for  the  leeward  fin  at  <t>£  =  45° 
are  shown  in  parts  (f)  of  the  three  figures  for  Mm  =  1.6, 

2.36,  and  3.7.  The  predictions  are  good  only  for  6  =  0  at 
low  angles  of  attack,  but  otherwise  are  not  satisfactory. 

Some  calculations  have  been  made  with  the  supersonic 
panel  method  of  reference  12  at  =  1.60  to  study  further 
the  problems  of  predicting  fin-normal  forces  and  hinge  moments 
for  the  subject  wing-body  combination.  In  these  calculations 
the  fins  were  treated  as  lifting  surfaces,  with  seven  chordwise 
and  five  spanwise  rows  of  panels  on  each  fin.  Calculations 
were  made  for  <J>f  =  45°  and  135°,  for  6  =  0°  and  -15°,  and 
for  a  =  0°,  15°,  and  20°,  a  total  of  12  cases.  The  results 
of  these  calculations  are  compared  with  the  data  of  reference 
39  in  Figure  37.  What  is  immediately  clear  is  that  for  nearly 
all  cases  the  fin  normal-force  coefficient  and  the  fin  bending- 
moment  coefficient  are  predicted  well  by  the  panel  theory  which 
uses  the  same  vortex  strengths  and  positions  as  the  engineering 
prediction  method.  However,  hinge-moment  coefficients  are  not 
satisfactorily  predicted.  The  reason  that  the  hinge-moment 
coefficients  are  not  well  predicted  is  that  the  panel  theory 


does  not  give  accurate  axial  center-of-pressure  positions. 

An  important  reason  for  this  shortcoming  is  that  it  does  not 
account  for  second-order  effects  of  thickness  on  x/cr. 

The  panel  method  yields  good  normal-force  coefficients 
at  =  45°  and  135°  even  at  angles  of  attack  up  to  20°  where 
the  engineering  prediction  method  does  not.  It  is  natural, 
then,  to  try  to  calculate  the  hinge  moments  by  a  hybrid  method 
using  the  panel-method  normal-force  coefficient  and  the 
experimental  values  of  x/cr  for  the  wing  alone.  This  is 
certainly  promising  for  the  windward  fins  when  no  vortex 
effects  of  any  significance  occur.  Consider  now  Figure  59 
which  is  a  reproduction  of  Figures  54(d),  (e) ,  and  (f) . 

Figure  59(a)  shows  that  the  hinge-moment  coefficient  at 
<J>£  =  135°  and  6  =  -15°  is  much  better  predicted  by  the  hybrid 
method  than  by  the  engineering  prediction  method.  Excellent 
agreement  is  still  not  achieved  at  large  angles  of  attack 
probably  because  of  the  unporting  effect  discussed  in  connec¬ 
tion  with  Figures  55  and  56. 

In  principal  the  hybrid  method  as  used  for  the  windward 
fins  does  not  apply  to  the  leeward  fin  (<j>f  =  45°)  because 
the  fin  normal-force  increment  due  to  body  vortex  effects 
acts  at  a  different  position  than  that  of  the  normal  force 
without  vortices  present.  For  6=0  the  engineering  method 
does  a  good  job  of  calculating  normal  force  and  the  panel 
method  does  not  at  high  a,  see  Figure  59(b).  The  fact  that 
the  panel  method  does  a  good  job  of  calculating  normal  force 
at  6  =  -15°  is  thus  an  anomoly.  Hence,  do  not  think  the  hybrid 
method  is  reliable  for  leeward  fins  near  the  vortex. 

In  summary,  the  engineering  method  appears  to  predict 
hinge  moments  for  =  90°  well  enough  to  give  trends  with  a, 

6  ,  and  M^.  For  the  windward  fins  at  $£  =  135°,  a  hybrid 
method  based  on  the  panel  method  and  the  engineering  method 
yields  improved  accuracy.  At  high  angles  of  attack  with 
negative  deflection  (6  =  -15°)  and  large  angles  of  attack, 
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15°-20°,  satisfactory  predictions  have  not  been  made  for 
leeward  fins.  It  is  to  be  remembered  that  the  data  base  has 
been  extrapolated  both  in  aspect  ratio  and  Mach  number. 
Additional  comparisons  are  made  in  the  next  section. 

6.3  Comparisons  for  Army  Generalized  Missile 

6.3.1  Preliminary  remarks.-  A  considerable  number  of 
comparisons  have  been  made  for  the  Army  Generalized  Missile, 
the  body  of  which  is  shown  in  Figure  50.  We  will  consider 
results  for  two  canard  fins  and  one  tail  fin  shown  in 
Figure  51  from  tests  described  in  reference  42.  Fin  C4 
was  mounted  at  the  mid  hinge-line  position  and  at  the  aft 

hinge-line  position  for  the  canards  as  shown  in  Figure  50. 

Tail  fin  T2  was  mounted  at  the  mid  hinge-line  tail  position. 
Comparisons  will  be  made  for  the  tail  fins  with  and  without 
the  canard  fins  present. 

It  is  helpful  to  tabulate  the  various  conditions  for 
which  comparisons  have  been  made  before  getting  into  the  details. 
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It  was  not  possible  to  deflect  the  tail  fins  so  that  is 
zero  in  all  cases. 

The  geometric  properties  of  the  fins  are: 

S„  c.  c  s„  -  a 
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Both  canard  fins  are  outside  the  aspect-ratio  boundaries  of 
the  fin-on-body  data  base  although  C4  requires  only  a  minor 
extrapolation . 

6.3.2  Wing-alone  normal-force  coefficients  and  axial 
center-of-pressure  positions.-  The  C4  fin  results  are  for 
of  1.5  and  2.0,  the  Cg  fin  results  for  0.8,  1.3,  and  1.75, 
and  the  T2  results  for  1.75.  The  data  base  for  wing-alone 
normal-force  coefficient  in  reference  1  was  used  to  predict 
CN  for  the  fins  in  the  following  steps  for  fin  C^. 

Step  1:  Linearly  extrapolate  in  1/7R  in  the  X  =  1/2 

data  to  get  CN  at  JR  =  2. 31,  X  =  1/2  for  each  a. 
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Step  2:  Linearly  interpolate  in  X  in  the  JR  = 1  data  to 
get  a  ACN  due  to  going  from  X=0.5  to  X  =  0.3, 


ACN  ”  (V-  0?5  )  {CnT14  ‘  CNt15) 


Step  3:  Add  ACN  to  Cn2  31  1/2  to  9et  CN  at  ®  =2.31 
and  X  =  0.3. 


Step  4:  Plot  Cm-  A  versus  1/M  for  M  =  1.2,  2,  and 
“2.31,0.3  00  00 

3  for  all  a's. 

Step  5:  Read  desired  values  of  CN  from  plot  for  Mro  =  1.5 
and  2.0 

Fin  Cg  is  in  the  wing-alone  and  fin-on-body  data  bases 
of  reference  1  for  M  =0.8  and  1.3.  However,  we  have  used 
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the  above  procedure  to  get  the  desired  wing  alone  normal-force 
coefficients  by  extrapolating  from  PR  =2  to  PR  =3.53.  The 
extrapolated  normal-force  curve  for  =  1.3  is  compared 
in  Figure  60  with  the  data  given  in  Table  3  of  reference  1 
to  give  an  idea  of  the  errors  arising  in  extrapolation  from 
PR  =2  to  PR  =3.53.  The  CN  curve  for  =  0.8  was  estimated 
from  the  data  of  reference  6,  but  because  the  nonlinearities 
at  transonic  speeds  are  wing-section  dependent,  the  curves  did 
not  compare  well. 

Fortunately,  fin  T2  is  within  the  wing-alone  data  base  in 
reference  1,  and  its  normal-force  curve  at  M  =1.75  was 
obtained  following  the  five  steps  of  the  procedure  outlined 
above . 

Now  with  respect  to  axial  center-of-pressure  positions 
for  the  wings  alone,  the  following  procedure  was  used  in  all 
cases . 

Step  1:  Determine  the  center-of-pressure  position  from 
DATCOM  based  on  a  lifting  surface;  (x/cr)LT. 

Step  2:  Determine  the  change  in  center-of-pressure  posi¬ 
tion  due  to  thickness  A(x/cr)t  for  angles  of 
attack  up  to  that  for  shock  detachment  by  using 
the  program  described  in  Appendix  B. 

Step  3:  Add  the  quantities  in  steps  1  and  2  to  get  the 
wing  alone  axial  center-of-pressure  position. 


Step  4:  Determine  the  axial  center-of-pressure  position 
at  a  =  45°  from  the  curves  of  Figures  34  or  35. 


Step  5:  Fair  a  curve  between  the  results  of  steps  3 
and  4 . 

The  (x/cr)w  curves  for  wing  C4  at  Mo  =  1.5  and  2.0  based  on 
this  method  are  shown  in  Figure  61;  those  for  wing  Cg  in 
Figure  62  for  Mach  numbers  of  1.3  and  1.75;  and  that  for  wing 
T2  at  =  1.75  in  Figure  63.  We  now  make  comparisons  between 
experiment  and  theory  in  the  order  listed  in  Section  6.3.1. 

6.3.3  Comparisons  for  fin  C4.-  The  comparisons  between 

experiment  and  prediction  for  normal  force  and  hinge  moment  for 

fin  C4  are  contained  in  Figure  64 .  The  comparisons  were  made 

for  M  =1.5  and  no  roll  with  the  horizontal  canards  undeflected 
00 

and  deflected  +15°.  The  notation  for  the  canard  fins  F^ ,  F 2, 
F3,  and  F4  follows  the  convention  below  looking  upstream. 


4>  =  0°  4>  =  45° 

We  note  that  the  fin  normal-force  coefficients  are  well 
predicted  by  the  method.  The  hinge-moment  coefficients  are 
small  and  the  trends  with  angle  of  attack  are  fairly  well 
predicted  even  though  the  data  and  predictions  have  opposite 
signs.  Inquiring  more  closely  into  the  agreement  between 
data  and  prediction,  we  note  that  the  normal- force  coefficients 
at  high  angles  of  attack  are  well  predicted  for  both  6=0° 


and  6  =  15°  so  that  any  disagreement  for  hinge-moment  coef¬ 
ficient  is  due  mainly  to  axial  center-of-pressure  position. 

For  6  =  15°  and  o  =  12.5°  the  difference  between  data  and 
prediction  represents  an  error  in  axial  center-of-pressure 
position  of  about  3  percent  of  the  root  chord  for  fin  F4  and 
about  5  percent  for  fin  F2«  It  is  noted  that  the  experimental 
difference  of  2  percent  between  the  two  fins  is  not  much  less 
than  the  difference  between  experiment  and  theory  for  fin  F4 . 
The  agreement  might  thus  range  from  good  to  fair. 

Results  are  shown  in  Figure  65  for  =  2.0  similar  to 
those  in  Figure  64  for  =  1.5  with  the  same  accuracy  of 
comparison.  There  are  no  significant  body  vortex  effects  for 
the  range  of  angles  of  attack  tested. 

Comparison  for  a  roll  angle  of  45°  and  Mach  number  of  1.5 
with  no  fin  deflection  is  shown  in  Figures  66  and  67.  It  is 
noted  that  data  are  available  on  the  fin  C4  quantities  with 
and  without  tail  fins  present,  and  the  duplicate  data  are 
shown  in  the  figures.  For  the  upper  fins  the  normal-force 
data  in  Figure  66 (a)  do  not  show  good  agreement  for  left  and 
right  fins.  For  the  lower  fins  all  the  data  agree  well  in 
Figure  66(b)  and  are  in  good  agreement  with  prediction.  The 
hinge-moment  results  of  Figure  67  are  also  satisfactory. 

Figures  68  and  69  give  normal-force  coefficient  and  hinge 
moment  results,  respectively,  for  =  1.5  and  <j>c  =  45°  with 
all  four  fins  positively  deflected  15°.  For  the  upper  two 
fins,  F ^  and  F2,  the  experimental  normal-force  coefficients 
and  the  hinge-moment  coefficients  do  not  agree  well  but  for 
the  lower  fins  the  data  are  very  repeatable.  The  comparisons 
between  experiment  and  theory  are  judged  to  be  good. 

In  Figures  70  and  71  normal- force  and  hinge-moment 
coefficients  are  presented  for  Mro  =  2.0,  ac = 45°  with  all  fins 
undeflected.  All  results  are  in  good  agreement  and  repeatable 
except  for  the  F.  and  F,  normal-force  coefficients.  Despite 


the  fact  that  and  F2  are  different,  they  repeat  tail  on 
and  tail  off.  However,  the  theory  is  well  below  the  data  for 
either  or  them.  This  tendency  was  also  seen  for  fin  Tg  at 
4>f  =  4  5°  in  Figure  57.  The  comparable  quantities  are: 


M 

OO 

*0 

JR 

X 

a/s 

m 

Fin 

T  • 

B' 

2.00 

45° 

2.812 

0.423 

0.25 

Fin 

c  • 

4  • 

2.36 

450 

2.31 

0.30 

0.54 

We  do  not  know  the  reason  for  this  discrepancy  but  suspect 
that  the  values  of  normal-force  coefficient  in  the  data  base 
for  any  of  the  three  fins  or  two  Mach  numbers  which  yield  the 
predicted  normal-force  coefficients  are  low.  The  data  in  the 
fin-on-body  data  base  was  taken  with  the  fin  trailing  edge 
aligned  with  the  base  of  the  body.  At  angle  of  attack  and 
=  45°  some  possible  base  effect  on  fin  normal  force  might 
occur.  Both  fins  Tfi  and  C4  had  at  least  some  afterbody  behind 
their  trailing  edges.  The  fins  in  the  tri-service  data  base 
are  being  tested  with  some  afterbody,  and  these  data  should 
resolve  this  question. 

Data  similar  to  those  for  Figures  70  and  71  are  shown  in 
Figures  72  and  73  except  will  all  fins  deflected  15°  rather 
than  0°.  Again  the  repeatability  of  the  data  and  the  agreement 
with  prediction  are  good  except  for  the  normal-force  coef¬ 
ficient  on  fins  F^  and  F2.  While  CNCl  and  CNC2  are  repeatable 
tail-on  and  tail-off,  they  are  not  equal  as  required  by  lateral 
symmetry.  If  F2  is  correct,  then  the  agreement  with  prediction 
is  fair. 

6.3.4  Comparisons  for  fin  C6--  As  mentioned  previously, 
the  experimental  fin-on-body  normal-force  coefficients  at 
M  =0.8  and  M  =1.3  are  included  in  the  data  base  used  in  the 
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predictive  method.  Nevertheless,  comparisons  are  made  between 
data  and  predictions  for  Cg  because  hinge  moments  are  not  in  the 


data  base  and  because  =  1.75  predictions  will  be  made. 

The  wing-alone  normal-force  curve  has  been  estimated  as 
described  previously. 

We  first  look  at  normal- force  and  hinge-moment  coef¬ 
ficients  in  Figure  74  for  =  1.3,  <J>C  =  0°,  and  6.^  =  6^  =  0 
and  15°.  The  equality  of  the  normal-force  coefficients  left 
to  right  is  good,  but  for  the  hinge-moment  coefficients  it 
could  be  better.  However,  the  predictions  for  normal-force 
coefficients  are  very  good  at  6  =  0°  and  6  =  15°,  and  the 
hinge-moment  predictions  match  the  data  well  in  the  mean. 

It  is  to  be  noted  that  the  significant  error  in  estimating 
the  wing-alone  normal-force  coefficient  shown  in  Figure  60 
apparently  has  little  effect  on  the  accuracy  of  the  predicted 
fin-on-body  normal  force  since  is  calculated  internal  to 
the  program  from  the  fin-on-body  data  base. 

Data  for  M  =  1.75  at  <J)  =0  with  6_  =  6.  =  0  and  15° 
are  exhibited  in  Figure  75.  The  agreement  between  the  data 
for  left  and  right  fins  and  the  comparison  of  experiment 
and  prediction  are  comparable  to  those  for  =  1.30  in  the 
previous  figure. 

Data  for  =  1.30  at  4>c  =  45°  with  6-^  =  62  =  ^3  =  ^4  =  0 
and  15°  are  given  in  Figure  76.  The  data  for  the  upper  fins, 
F-^  and  Fjf  show  some  difference  left  to  right  at  6  =  0  and 
6  =  15°  for  both  normal-force  and  hinge-moment  coefficients. 
However,  the  agreement  between  experiment  and  theory  is  con¬ 
sidered  fairly  good  up  to  a  =  20°.  The  lower  fins,  Fj  and  F4 
show  good  repeatability  of  the  data  left  to  right  for  both 
6=0  and  6  =  15°  for  both  coefficients.  The  comparison 
between  experiment  and  theory  are  good  for  normal- force  coef¬ 
ficient  and  for  hinge  moment  for  6=0.  For  6  =  15°  the 
hinge-moment  agreement  is  only  fair.  It  is  suspected  that 
body  vortex  effects  are  not  important  in  this  case,  and 
unporting  may  be  a  contributing  factor  through  its  effect  on 
axial  center-of-pressure  position. 


Data  for  Mw  =  1.75  are  given  in  Figure  77  paralleling 
that  for  =  1.3  in  Figure  76.  The  agreement  of  the  normal- 
force  and  hinge-moment  coefficients  left  to  right  for  the 
upper  fins  F^  and  F2  is  not  good  for  6=0  nor  6  =  15°.  For 
6=0°  the  normal-force  prediction  for  F^  and  F2  is  low, 
but  at  6  =  15°  it  is  good  up  to  a  =  20°.  The  hinge-moment 
comparisons  for  6  =  0  is  fair  up  to  20°  and  for  6  =  15°  is 
good  in  the  mean.  For  the  windward  fins  F^  and  F^  the  data 
are  very  close  left  to  right  and  the  comparisons  between 
experiment  and  theory  are  good  for  all  cases. 

An  attempt  was  made  to  predict  the  hinge  moments  of 
fin  C6  at  Mot  =  0.8,  4>c  =  0,  with  6  =  0  and  6  =  15°.  However, 
the  attempt  was  not  successful.  The  wing-alone  characteristics 
for  Cg  had  to  be  estimated  from  the  Stallings-Lamb  data  base 
since  one  of  the  required  wings  alone  was  not  present  in  the 
Fidler-Baker  data  base.  However,  the  differences  in  wing 
sections  between  Cg  and  the  equivalent  Stallings-Lamb  wing 
were  too  large  to  yield  accurate  center-of-pressure  predictions 
at  M  =0.8. 

CO 

6.3.5  Comparisons  for  tail  fin  T2--  Tail  fins  T2 
exhibited  phenomena  different  from  canard  fins  C4  and  Cg 
because  they  are  mounted  well  back  on  the  body  where  vortex 
effects  are  significant.  Also  they  were  tested  with  and 
without  canard  fins  present  so  that  the  effects  of  wing-tail 
interference  on  hinge  moments  could  be  investigated.  The 
following  listing  shows  the  various  conditions  for  which  com¬ 
parisons  between  experiment  and  prediction  are  made.  The 
model  did  not  permit  tail  fin  deflection.  When  the  canard 
fins  were  deflected,  all  four  fins  were  positivelydef lected 
15°  at  <J>c  =  45°,  but  only  the  horizontal  canards  were  deflected 
at  <Pc  =  0°  for  the  cases  presented. 
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For  the  tail  fins  the  numbering  of  the  fin  has  been 
done  clockwise  rather  than  counterclockwise  as  for  the  canard 
fins.  Thus,  the  fin  numbering  for  the  4>  =0°  and  4>c  =  45° 
cases  is  as  follows. 


The  agreement  of  the  data  between  fins  F2  and  F^  (body- 
tail  configuration)  in  Figure  78  is  very  good,  and  the  agree¬ 
ment  between  data  and  theory  for  normal-force  coefficient 
is  also  very  good.  The  hinge-moment  coefficient  prediction  is 
only  fair  at  angles  of  attack  approaching  20°.  At  the  highest 
angle  of  attack,  the  error  in  center-of-pressure  position 
is  .05. 

For  4>c  =45°  and  no  canards  present.  Figure  79  also  shows 
excellent  agreement  between  the  measured  results  for  the 
two  upper  fins  and  for  the  two  lower  fins.  The  lower  fins, 

2  and  3,  show  very  good  agreement  between  experiment  and 
theory  for  normal  force,  and  the  upper  fins,  1  and  4,  show 
good  agreement  when  allowance  is  made  for  the  large  reduction 


in  normal  force  associated  with  the  body  vortices.  The  hinge- 
moment  coefficients  for  the  windward  fins  are  underpredicted 
in  about  the  same  way  as  for  the  horizontal  fins  in  Figure  78. 
This  result  is  not  a  vortical  effect.  The  leeward  fins 
show  a  substantial  reduction  in  hinge-moment  coefficients 
below  that  for  the  windward  fins  for  which  the  prediction  is 
fair. 

Figure  80  shows  the  normal-force  coefficient  and  hinge- 
moment  coefficients  when  canard  fins  are  added  to  the  body- 
tail  combination  which  is  treated  in  Figure  78.  Some  reduction 
in  normal  force  occurs,  but  not  as  much  as  predicted  at  angles 
of  attack  of  20°  or  greater.  The  measured  hinge  moments 
increased  only  very  slightly,  but  the  prediction  is  now  in 
better  agreement  with  the  data.  Vortices  from  the  canard  fins 
have  trailed  back  and  passed  over  the  tail  fins,  but  are 
fairly  far  above  them  before  the  vortex  strength  is  large. 
Therefore,  vortex  effects  are  small. 

In  Figure  81  the  same  results  are  shown  as  for  Figure  80 
except  that  the  forward  canards  have  now  been  deflected  posi¬ 
tively  by  15°.  Canard  vortices  thus  pass  below  the  horizontal 
tail  fins  near  zero  angle  of  attack.  However,  they  intersect 
the  tail  fins  at  an  angle  of  attack  near  5°  causing  the  bumps  in 
predicted  normal-force  and  hinge-moment  coefficient  curves.  At 
a  =  0  the  download  on  the  tail  fins  is  due  entirely  to  the  canard 
vortices  and  is  fairly  well  predicted.  However,  the  hinge-moment 
coefficient  predicted  at  a  =  0  has  the  wrong  sign.  This  result 
warrants  refinement  in  the  method  of  determining  vortex  effects 
on  hinge  moments.  Overall  the  trends  with  angle  of  attack  are 
fairly  well  predicted. 

In  Figure  82  data  and  theory  are  compared  for  tail-fin 
hinge  moments  with  the  model  rotated  to  the  <j>  =  45°  position. 

In  this  case  there  is  a  pair  of  canard  vortices  originating  on 
the  bottom  fins  which  does  not  intersect  the  lower  pair  of 
tail  fins  but  which  intersects  the  upper  pair  at  an  angle  of 
attack  near  a  =  13°.  The  bumps  in  the  predictions  for 


normal-force  and  hinge-moment  coefficients  are  a  manifesta¬ 
tion  of  this  fact  together  with  the  way  in  which  the  computer 
program  treats  vortices  in  close  proximity  to  the  fin.  The 
agreement  between  experiment  and  theory  for  the  upper  fins 
is  not  good.  The  lower  fins  are  not  intersected  by  canard 
vortices,  and  the  agreement  between  experiment  and  theory  is 
good  to  fair. 

In  Figure  83,  the  same  case  is  handled  as  for  Figure  82 
except  that  the  canard  fins  are  all  positively  deflected  15°. 
Strong  vortices  exist  even  at  a  =  0°.  The  agreement  between 
experiment  and  theory  in  this  case  is  unsatisfactory. 

To  summarize,  the  prediction  method  gives  fair  to  good 

results  for  the  body-tail  combination  at  <j>  =0  and  <f>  =45°. 

c  c 

For  <b  =0  and  addition  of  the  canards  with  or  without  deflec- 

c 

tion  of  the  horizontal  canards,  the  predictions  of  the  method 
are  fair  to  good.  For  =  45°  and  addition  of  the  canards 
with  6=0  the  predictions  are  not  good  and  deflection  of  all 
four  canards  makes  prediction  even  worse. 

6.3.6  Comparison  between  data  and  theory  for  effects  of 
roll  angle  on  hinge  moment.-  Sufficient  data  are  available  to 
see  how  the  method  predicts  normal- force  and  hinge-moment 
coefficient  for  a  roll  angle  range  of  180°  for  fixed  angles 
of  attack  for  both  canard  fins  and  tail  fins.  Since  all  four 
fins  are  instrumented,  such  data  can  be  generated  by  rolling 
the  model  only  45°.  Comparisons  have  been  made  for  the 
following  cases: 
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In  these  plots  <J>  =  0 

is  the  horizontal  fin,  <J>  =  90°  the  windward  fin,  and  <J>  =  -90° 

the  leeward  fin.  The  data  on  hinge  moments  of  two  different 
fins  at  the  same  condition  give  a  measure  of  the  minimum  error 
of  the  data. 

Figure  86  for  M^  =  1.75  shows  an  underprediction  of  the 
leeward  fin  normal-force  coefficient  for  both  10°  and  20°. 

The  prediction  involves  extrapolation  both  in  Mach  number 
(from  1.3)  and  aspect  ratio.  Figure  87  exhibits  good 
comparison  between  data  and  theory  for  hinge-moment  coeffi¬ 
cient  at  a  =  10°,  but  underprediction  of  the  data  for  <j>c  <  -30° 
for  a  =  20°.  The  discrepancies  here  are  not  vortex  effects 
since  they  are  small  for  the  canard  fins. 

For  tail  fin  T2  the  normal-force  coefficients  in  Figure  88 
are  well  predicted  including  the  severe  body  vortex  effects 
exhibited  at  a  =  20°.  Despite  good  normal-force  prediction, 
the  hinge  moment  is  consistently  underpredicted  at  both  angles 
of  attack  in  Figure  89. 

7 .  CONCLUDING  REMARKS 

A  comprehensive  engineering  hinge-moment  predictive 
method  for  cruciform  missiles  with  all-movable  controls  has 
been  formulated  in  preliminary  form,  and  a  series  of  compari¬ 
sons  have  been  made  between  experiment  and  prediction  to 


determine  the  range  of  applicability  of  the  method  together 
with  its  shortcomings.  In  addition,  the  applicability  of 
a  supersonic  panel  method  for  cruciform  missiles  has  been 
investigated  with  respect  to  calculation  of  hinge  moments 
for  all-movable  controls.  The  engineering  method  in  principle 
is  applicable  from  transonic  to  high  supersonic  speed  for 
wing-body  or  wing-body-tail  combinations  employing  cruciform 
all-movable  controls.  The  method  applies  to  large  angles  of 
attack  and  control  deflections,  the  sum  of  which  can  approach 
50°.  The  panel  method  is  applicable  in  the  supersonic  range 
but  covers  a  much  broader  spectrum  of  configurations  than  the 
engineering  method.  It  is,  however,  limited  to  about  20°  of 
angle  of  attack  plus  control  deflection. 

For  accurate  hinge-moment  prediction,  it  has  been  found 
that  accurate  knowledge  is  required  of  the  fin  axial  center- 
of-pressure  positions  —within  about  2  percent  of  the  root  chord. 
Large  changes  in  axial  center-of-pressure  position  occur  with 
angle  of  attack  and  Mach  number,  but  existing  analytical  methods 
are  not  able  to  predict  them  within  the  desired  accuracy.  A 
hybrid  method  of  determining  axial  center-of-pressure  position 
of  wings  alone  has  been  developed  for  wings  having  attached 
leading-edge  shock  waves  when  viewed  in  planes  normal  to  the 
wing  planform  and  parallel  to  the  root  chord.  This  method 
provides  essential  center-of-pressure  information  which  is  used 
in  both  the  engineering  method  and  the  panel  method  for  estimat¬ 
ing  hinge  moments.  In  the  low  supersonic  region  (M^  <  1.2)  and 
the  high  subsonic  region  (Mro  >  0.8),  the  fin  center-of-pressure 
position  is  a  function  of  wing  airfoil  section  in  a  way  which 
cannot  be  currently  estimated.  Accordingly,  in  these  regimes 
existing  data-base  information  for  wing-alone  center-of-pressure 
position  must  be  used.  This  fact  limits  the  configuration 
range  that  can  be  treated. 

Comparisons  have  been  made  for  a  number  of  configurations 
over  a  range  of  operating  parameters  to  determine  the 
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limitations  of  the  engineering  predictive  method.  Consider 
first  the  comparisons  for  a  cruciform  body-tail  combination 
with  fins  Tfi  having  A*  =2.812  and  X  =  0.423.  Comparisons 
were  made  for  the  following  ranges  of  parameters: 

M  :  1.6,  2.35,  3.7 

00 

a:  0°  £  a  £  20° 

6:  0°,  -15° 

<J)f :  45°,  90°,  135° 

For  the  unrolled  case  (4>^  =  90°)  the  prediction  and  data  com¬ 
pared  well  for  the  ranges  of  the  first  three  parameters.  With 
a  pair  of  fins  on  the  windward  side  ( ^  =  135°)  it  was  found 
that  some  times  the  engineering  method  predicted  the  fin  normal 
force  inaccurately,  but  the  combined  use  of  the  panel  and 
engineering  method  improved  accuracy  of  prediction.  For  the 
leeward  fins  (4>f  =  45°)  ,  the  fin  hinge-moment  predictions  were 
usually  poor.  In  these  cases  strong  body  vortices  appeared 
at  the  high  angles  of  attack. 

A  second  set  of  systematic  comparisons  was  made  for  a 
body  with  canard  fins  (C^)  having  an  aspect  ratio  of  2.36  and 
taper  ratio  of  0.30.  The  comparisons  were  made  for  the 
following  ranges  of  parameters: 

M  :  1.5,  2.0 

OO 

a:  0°  £  a  £  20° 

6:  0°,  +15° 

4>f:  45°,  90°,  135° 

In  some  cases  the  left  and  right  firs  did  not  give  equal 
measured  results  when  they  should  have  by  symmetry.  However, 
one  set  of  results  agreed  fairly  well  with  prediction  in  all 
cases  but  one.  In  this  one  case  of  approximately  40  compari¬ 
sons,  significant  disagreement  between  prediction  and  data 

occurred.  This  is  a  case  for  M  =  2.0  at  45°  roll  with  6=0 

00 
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where  the  normal-force  coefficient  for  the  leeward  fins  is 
significantly  underpredicted.  Figure  70(a).  A  conjecture  con¬ 
cerning  the  cause  is  given  in  the  text.  Generally,  the  good 
agreement  for  this  case  between  prediction  and  experiment 
can  be  ascribed  to  the  absence  of  any  significant  body  vortex 
effects  in  contrast,  for  example,  to  the  previous  case. 

The  third  set  of  comparisons  was  made  for  another  body- 
canard  fin  combination  with  fins  (Cg)  of  aspect  ratio  of  3.53 
with  X  =  .06.  The  ranges  of  parameters  are: 

M  :  1.3,  1.75 

00 

a:  0°  <  a  <  20° 

6:  0°,  15° 

<J)f :  45°,  90°,  135° 

Essentially  the  same  results  were  found  for  the  previous  body- 

canard  combination.  An  attempt  to  predict  hinge  moments  at 

M  =0.8  was  unsuccessful  because  of  the  lack  of  a  reliable 
00 

way  of  predicting  fin  or  wing-alone  axial  center-of-pressure 
position  at  subsonic  speeds. 

The  fourth  set  of  comparisons  was  made  for  a  body-tail 
combination  employing  tail  fins  T2  to  which  canard  fins  C4 
were  then  added  to  see  how  well  the  effect  of  wing-tail  inter¬ 
ference  on  hinge  moments  could  be  predicted.  The  tail  fins 
have  an  aspect  ratio  of  1.333  and  a  taper  ratio  of  0.50.  The 
parameters  are: 

M  :  1.75 

00 

a:  0°  £  a  <.  20° 

6c:  0°,  15°  (canards  only) 

4>f:  45°,  90°,  135°  (tail  fins  inline) 

First  the  comparisons  between  prediction  and  data  for  the 
body-tail  combination  is  discussed  as  a  basis  for  evaluating 
the  subsequent  effects  of  adding  canard  fins.  For  =  90° 
and  135°  the  fin  normal  forces  are  predicted  quite  closely. 


The  leeward  fin  normal-force  coefficient  is  also  predicted 
well  considering  the  fact  that  at  a  =  20°  it  is  less  than  30 
percent  of  that  for  the  windward  fins.  The  hinge-moment 
coefficients  are  underpredicted  by  about  25  percent  as  a 
rough  average. 

Adding  the  undeflected  canard  fins  at  =  1.75  does  not 
have  a  large  effect  on  the  tail  fins,  and  the  agreement  between 
prediction  and  data  for  the  hinge  moment  is  improved.  Deflect¬ 
ing  the  horizontal  canards  15°  positive  cause  nonlinearities 
in  normal  force  and  hinge  moment  due  to  intersection  of  the 
tail  fin  by  the  canard  vortices,  nonlinearities  which  do  not 
appear  in  the  data.  However,  the  agreement  between  prediction 
and  data  are  still  fair.  When  the  configuration  is  rolled 
45°,  the  effects  of  the  canards  on  the  tail  fin  hinge  moments 
are  not  well  predicted,  with  or  without  canard  deflection. 

The  last  set  of  comparisons  were  made  to  see  how  the 
effects  of  roll  angle  on  normal  force  and  hinge  moment  were 
predicted  for  a  roll  angle  range  of  180°.  The  first  compari¬ 
sons  for  the  body-canard  (Cg)  configuration  for  =  1.3 
exhibited  good  agreement  between  experiment  and  theory  at 
both  a  =  10°  and  20°  over  the  roll  angle  range.  At  =  1.75 
there  was  a  tendency  to  underpredict  the  normal  force  on  the 
leeward  fins  at  both  a  =  10°  and  20°.  This  tendency  was 
reflected  in  the  hinge-moment  coefficient  at  a  =  20°  only. 

To  see  the  effect  of  body  vortices  on  the  preceding  compari¬ 
sons,  similar  comparisons  were  made  for  the  body-tail  (T2) 
configuration.  The  normal-force  comparisons  are  quite  good, 
but  the  hinge-moment  coefficients  are  underpredicted. 

These  comparisons  indicate  the  following  conclusions: 

1.  The  method  is  applicable  to  the  range  1.3  <  <  3.7. 

2.  It  is  not  applicable  with  accuracy  to  leeward  fins 
near  strong  body  vortices. 


V, 


3.  A  hybrid  theory  is  applicable  to  windward  fins. 


4.  It  will  handle  wing-tail  interference  effects  for 
inline  canards  and  tails  with  no  roll  angle  both 
with  and  without  canard  deflection. 

8 .  RECOMMENDATIONS 

As  a  result  of  the  present  work,  several  problem  areas 
have  arisen  in  which  further  work  should  result  either  in 
increasing  the  accuracy  of  hinge-moment  estimation  or  in 
extending  the  range  of  applicability  of  the  present  method. 

1.  For  Mach  numbers  less  than  the  Mach  number  of  shock 
detachment  for  zero  degrees  angle  of  attack  for  the  wing-alone, 
data  have  exhibited  significant  effects  of  wing  thickness  on 
axial  center-of-pressure  location.  No  method  is  available  for 
calculating  such  effects.  If  one  were  available,  the  range  of 
the  hinge-moment  prediction  method  could  be  significantly 
increased. 

2.  When  strong  body  or  canard  vortices  pass  in  close 
proximity  to  tail  fins,  the  present  calculative  model  for 
hinge  moments  seems  to  be  inaccurate.  This  could  result  from 
inadequacy  of  the  vortex  itself  or  in  the  method  of  determining 
its  effect  on  the  fin.  Resolution  of  this  problem  may  depend  on 
a  good  experiment  measuring  both  vortex  properties  and  fin 
pressure  distributions.  While  this  has  been  done  at  subsonic 
speeds,  a  parallel  program  is  needed  at  supersonic  speeds. 

Then  an  improved  analytical  method  should  be  developed. 

3.  For  wings  of  taper  ratio  0.5,  the  strip  theory/shock- 
expansion  method  of  determining  the  effects  of  thickness  on 
axial  center-of-pressure  position  was  not  accurate  at  low 
supersonic  Mach  number  even  though  the  leading-edge  shock  was 
attached.  It  nevertheless  was  accurate  for  taper  ratios  of 

0  and  1.0.  Further  study  of  the  problem  is  required  to  resolve 
this  apparent  anomoly. 
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4.  When  the  tri-service  data  base  has  been  measured,  the 
present  program  based  on  MISSILE  2  should  be  expanded  to 
encompass  the  greater  ranges  of  aspect  ratio  and  Mach  number. 

Also  these  data  will  yield  systematic  values  of  fin  effectiveness 
as  a  function  of  the  various  parameters  as  well  as  fin  inter¬ 
ference  factors,  This  information  should  be  used  to 

improve  the  accuracy  and  range  of  applicability  of  the  present 
program. 

5.  Systematic  data  are  needed  on  all-movable  controls  of 
general  planform,  particularly  ones  with  swept  trailing  edges 
so  that  the  method  outlined  in  this  report  to  calculate  their 
hinge  moments  can  be  properly  tested. 

6.  In  several  instances  it  has  been  observed  that  zero 
fin  normal  force  and  zero  fin  hinge  moments  did  not  occur  at 
the  same  conditions.  This  fact  implies  that  some  phenomena 
exist  which  exhibit  induced  camber  effects.  A  special  study 
to  understand  and  predict  these  effects  seems  worthwhile. 

7.  There  is  evidence  to  suggest  that  the  length  of  the 
afterbody  behind  the  trailing  edge  of  an  all-movable  control 
can  have  a  significant  effect  on  its  normal  force  and  hinge 
moment  in  certain  circumstances.  An  experimental  investiga¬ 
tion  of  this  question  is  recommended. 
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Figure  1.-  Summary  of  characteristics  of  wings 
of  Stallings  -  Lamb  data  base. 
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Effect  of  v/edge  angle  on  center-of -pressure  position  (positive  forward 
from  the  50  percent  chord  location)  for  three  airfoils  of 
the  same  thickness  ratio;  M_  =  4.50. 
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Figure  4.  -  Effect  of  Mach  number  on  center-of-pressure  position  (positive  forward 

from  50  percent  chord  location);  airfoil  2. 
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Figure  15.-  Comparison  of  axial  center-of-pressure  positions  for 
wing-alone  and  fin-on-body  at  6=0°;  PR  =  1.0,  X  =  0. 
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Figure  17.-  Continued 
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Figure  19.—  Comparison  of  wing— alone  nonlinear  center-of— pressure  travel 
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(b)  a  -  20°,  M  =  2. 1C,  y/(b/2)  =  0 


Figure  23.-  Longitudinal  pressure  distributions 
on  an  aspect  ratio  2.0  rectangular  wing. 
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Figure  28.-  Longitudinal  center-of-pressure  locations  of  an 
aspect  ratio  1.0,  taper  ratio  0.5  cropped  delta  wing. 
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Figure  30.-  Longitudinal  center-of-pressure  locations  of  an 
aspect  ratio  4.0,  taper  ratio  0.5  cropped  delta  wing. 


Figure  30.-  Concluded 
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Figure  31.-  Longitudinal  center-of-pressure  locations  of  an 
aspect  ratio  0.5  rectangular  wing. 
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Figure  32.- 


Longitudinal  center-of-pressure  locations  of  an 
aspect  ratio  1.0  rectangular  wing. 
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Figure  38.-  Application  of  panel  method  to  determine  axial  center-of-pressure 
position  of  fin  T„  in  presence  of  the  body. 
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Figure  42.-  Accuracy  of  linear  extrapolation  of  wing-alone 
axial  center-of-pressure  position  to  Mach  number 
greater  than  three. 
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Figure  45.-  Data  illustrating  accuracy  of  predicting  wing 
axial  center-of-pressure  position  for  JR  >  2  from 
data  for  PR  =  1  and  JR  =  2. 


axial  center-of-pressure  position,  x/c 


normal-force  coefficient 


Data, 

ref.  3 

a 

o 

o 

0 

JL 

□ 

20° 

0 

OJ 

o 

0 

A 

_ 1 _ 

40° 

- 1 _ 1 _ 

X,  taper  ratio 

(a)  M  =1.6;  iR  =  2 

00 

Figure  46.-  Variation  of  wing  normal-force  coefficient  with 
taper  ratio  for  iR  =  2  wings  of  Stalling-Lamb  data  base. 
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Figure  46.-  Concluded. 
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Figure  49.-  Details  of  model.  All  linear  dimensions  are  in  centimeters 


Note  1:  All  dimensions  are  in  inches. 

Note  2:  Leading-  and  trailing- edge  radii  arc  0.020. 
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Figure  51.-  Continued 


Figure  53.-  Axial  center-of-pre»sure  position  for  wing 
alone  TB  at  several  Mach  numbers  by  interpolating  in 
Stallings-Lambdata  base. 
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Figure  54.-  Continued 
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Figure  54.-  Continued. 
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Figure  57.-  Continued. 
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Figure  58.-  Comparison  of  prediction  and  measurement  for  fin 


TB  characteristics  at  =  3.70;  Lamb-Trescot  data. 
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Figure  58.-  Continued 
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Figure  58.-  Continued 
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Figure  59.-  Comparison  of  engineering  prediction  method  and 
hybrid  method  with  measurments  for  fin  TR  characteristics 
at  ■  1.6;  Lamb-Trescot  data. 
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Figure  59.-  Continued. 
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Figure  63.-  Predicted  center-of-pressure  location  for  tail 
Tj  alone;  JR  =  1.33,  X  =  0.5. 
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Figure  72.-  Canard  C*  fina  normal-force  coefficients; 
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Figure  75.-  Concluded. 
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Figure  79.-  Characteristics  of  tail  T2  in  the  presence 
of  the  body;  M»  *  1.75,  4C  ■  45° 
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Figure  81-  Characteristics  of  tail  T2  in  presence  of  canard 
C6;  *  1.75,  4>c  =  0*,  »  S3  =  S2  =  &A  =  15°. 
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(b)  Hinge-moment  coefficient 
Figure  82.-  Concluded 
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Figure  84.-  Effect  of  roll  angle  on  normal-force  coefficient  for 
canard  C-  mounted  on  body  BN,;  M  =  1.3,  6.=6  =6  =6  .  =  0° 
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APPENDIX  B 


EFFECT  OF  THICKNESS  ON  AIRFOIL  CENTER-OF-PRESSURE 
POSITION  IN  ATTACHED  SHOCK  REGION 

B.l  CALCULATION  OF  PRESSURE  DISTRIBUTION  ON  AN 
AIRFOIL  USING  SHOCK-EXPANSION  THEORY 

The  method  of  calculating  the  pressure  distribution  on  a 
symmetric  airfoil  using  shock-expansion  theory  will  be  described. 
This  method  is  valid  up  to  the  angle  of  attack  at  which  the 
leading-edge  shock  wave  becomes  detached. 

Figure  Bl  shows  a  symmetric  airfoil  section  whose  surface 
is  made  up  of  a  series  of  straight  line  segments.  The  value 
of  x,  xfi,  at  which  the  segment  ends  is  specified  as  is  the 
surface  angle  0n.  On  the  front  of  the  airfoil  0n  is  positive. 
Downstream  of  the  point  of  maximum  thickness  0n  is  negative. 

On  the  lower  surface,  the  flow  in  each  region  is  characterized 
by 

Cpo  =  pressure  coefficient,  (p /  -  p  )/q 

cn  •'-n  00  00 

Mp  =  Mach  number 
-<-n 

=  local  Prandtl-Meyer  angle 

Similarly,  the  flow  in  each  region  on  the  upper  surface  is 
characterized  by 

Cp  =  pressure  coefficient,  (pu  -  p  )/q 

un  00  00 

Mn  =  Mach  number 
un 

vUn  =  local  Prandtl-Meyer  angle 

The  calculation  of  the  N  values  of  Cp^  and  the  N  values  of  Cp 
will  now  be  described.  Most  of  the  equations  used  are  taken 
from  reference  Bl. 


B.1.1  Lower  Surface  of  Airfoil 


On  the  lower  surface  in  region  1  the  free-stream  flow  is 
deflected  through  an  angle  6  where 


6  «  o  +  6X  (B-l) 

To  calculate  the  flow  quantities  in  region  1,  the  shock  wave 
angle,  6S,  which  is  a  function  of  and  6  must  be  determined. 
This  angle  can  be  determined  by  an  iterative  solution  of  the 
following  equation  (y  =  1.4  for  air) 


ctn6 


tan0s 


<Y  +  DM* 

2 (M?sin20  -  1) 

oo  S 


1 


(B-2) 


This  equation  is  double  valued  in  0  .  Also,  for  a  given  value 

s 

of  M  ,  there  is  a  value  of  6  above  which  no  solution  can  be 
00 

found,  the  case  of  a  detached  shock  wave.  The  maximum  value 
of  the  wedge  angle,  6,  which  will  allow  an  attached  shock  occurs 
when,  from  Equation  (B-2) , 


xrs  -  0  <B-3> 

Therefore,  the  maximum  value  of  0g  for  an  attached  shock  wave 
for  a  given  Ma  can  be  found  by  differentiating  Equation  (B-2) 
and  setting  the  derivative,  d6/d0  ,  to  zero.  If  this  is  done 


sin0 


s_ 


max 


-[4-(y+1)M2]  +  /(y+1)  V+8(y+1)  (y-1)M^  +  16(y+1) 


11/2 


4yM^ 


(B-4 ) 


Using  0q  and  M  ,  Equation  (B-2)  can  be  used  to  find  the 
&max  " 


maximum  value  of  6  for  an  attached  shock  wave . 

The  procedure  to  be  used  in  solving  Equation  (B-2)  for  6 


for  given  values  of  M  and  6  is: 

oo 

1.  Determine  0Smax  usi-n9  Equation  (B-4) 


2.  Using  0S  in  Equation  (B-2)  compute  <5 
max 


max 


3.  If  6  <  6  solve  Equation  (B-2)  for  the  value  of  6 


max 

which  is  less  than  0 


Smax 


Pressure  coefficient: 


Pl1  “  poo  4(M2sin20  -  1) 

1  00  s 


(Y  +  1)M" 


(B-5) 


Mach  number: 


Mi 


(y+1) 2M^sin20s  -  4 (M2sin28s  -  1)  (yM2sin26s  +  1) 


[2YM^sin2es  -  (y-1) ]  [ (Y-l)M2sin20s  +  2] 
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(B-6) 


Prandtl-Meyer  angle  for  : 


■y 


y:: 


■ .  i 


V 


.V< 


m 


>jS 


si 


V- 


Region  n: 

The  procedure  used  to  calculate  the  pressure  coefficients 
in  regions  2,  3,  . ...N  on  the  lower  surface  is  identical  so 
that  the  following  equations  are  written  for  region  n  where 
n  =  2,  3,  ... .N. 


The  calculation  is  repeated  sequentially  for  all  these  regions 

The  flow,  in  going  from  region  n-1  to  region  n,  expands 
through  angle  (0  ^  “  9n) •  Therefore,  the  Prandtl-Meyer  angle 

in  region  n  is 


=  v. 


n 


'n-1 


<9n-l  -  V 


(B-! 


With  this  angle  known,  the  Mach  number  in  region  n  can  be 
calculated  using  equations  given  in  reference  B2. 


M 


-2  -3 

1  +  1.3604  v  +  0.0962  v  -  0.5127  v 
-  I  zj 

n  1  -  0.6722  v  -  0.3278  v 


(B-l 


where 


The  ratio  of  static  pressure  to  total  pressure  in  region  n 
is  then 


pt,  (  i  s 

n 


and  the  pressure  coefficient  is 


(B-12) 


(B-13) 


since  pt  =  pt  .  Eauation  (B-8)  gives  (pt  /pt  )  and 

In  -C 1  “  L-C 1  oo 


Y_1 

00 


^oo  y  2  ^°° 
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A  11  _ 

2  oo  p 
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(B-15) 


After  the  calculations  given  by  Equations  (B-9)  through  (B-15) 
have  been  done  sequentially  for  n  =  2,3,...N,  the  complete 
pressure  distribution  on  the  lower  surface  has  been  determined. 


B.1.2  Upper  Surface  of  Airfoil 


Region  1: 


Two  separate  methods  are  used  in  calculating  the  pressure 
coefficient,  CPU^»  in  region  1  on  the  upper  surface  depending 
on  whether  the  angle  of  attack,  a,  is  less  than  or  greater  than 
the  surface  angle,  9^,  shown  in  Figure  Bl.  The  two  methods 
will  now  be  described. 


B-5 


a  <  e. 


Shock  Wave 


If  o  <  0^  the  free-stream  flow  is  deflected  through  an 


angle  where 


6  =  -(a  -  ex) 


(B-16) 


and  a  shock  wave  exists  on  the  upper  surface.  For  this  case  the 

procedure  described  for  region  1  on  the  lower  surface  is 

followed.  The  value  of  S  given  by  Equation  (B-16)  is  used  along 

with  the  free-stream  Mach  number,  M^,  in  Equation  (B-2)  to  find 

the  shock-wave  angle,  0_.  Then 

s 


Pressure  coefficient: 


\  -  Poo 


4  (M2sin20.,  -  1) 

co  S 


(y+dm; 


Mach  number: 


(y+1) 2M*sin20  -  4 (M2sin20  -  1) (yM2sin20s  +  1) 


[2yM~sin  0g  -  (y-1) ] [ (y-l)M;Sin“0g  +  2] 


Prandtl-Meyer  angle  for  Mu. : 


( y+1) M2sin20 


(y-1) M2sin20_  +  2 


2yM2sin20  -  (y-1) 

oo  5 


(B-20) 


a  >  0- 


If  a  >  0^  the  flow  expands  through  angle  (a  -  0^)  in  going 
from  the  free  stream  to  region  1  on  the  upper  surface.  There¬ 
fore,  the  Prandtl-Meyer  angle  in  region  1  is 


vux  =  v»  +  "  91> 


(B-21] 


where 


v  =,  pS.  tan’1,  E 
00  Vy-1  V  Y+- 


(M  -  1 )  -tan  /M 


-1  / m2  - 


(B-22! 


and  the  Mach  number  is 


1  +  1 . 3604  v  +  0.0962  v2  -  0.5127  v3 


(B-23 ; 


1  -  0.6722  v  -  0.3278  v‘ 


where 


(B-24) 


The  ratio  of  static  pressure  to  total  pressure  in  region  1  is 


1  + 


2zi 

2 


(B-25) 


and  the  ratio  of  total  pressure  in  region  1  to  free-stream 
total  pressure  is,  since  there  is  no  shock  wave. 


u. 


=  1.0 


(B-26) 


The  pressure  coefficient  in  region  1  on  the  upper  surface  is 
given  by 
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Ptul 

Poo 

Pt 

Pt 
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*01 
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ulJ 

00 

00 
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^00 


(B-27) 


The  quantities  P0()/Pt  and  qoo/pt  are  obtained  from  Equations 


(B-14)  and  (B-15) . 


Region  n: 

The  procedure  for  calculating  the  pressure  coefficients  in 
regions  2,  3,...N  on  the  upper  surface  is  the  same  as  that 
described  previously  for  the  same  regions  on  the  lower  surface. 
The  following  equations  are  written  for  region  n  where 
n  =  2,  3, . . ,N. 


The  flow  in  going  from  region  n-1  to  region  n  expands 

through  the  angle  (0  ,  -  6  ) .  The  Prandtl-Meyer  angle 

n-l  n 


in  region  n  is 


v  =  v  +  (0  ,  -  0  ) 

u  u  .  n-l  n 

n  n- 1 


(B-28 ) 


If  n=2,  either  Equation  (B-19)  or  (B-21)  is  used  to  determine 
vun_i  depending  on  whether  a  is  less  than  or  greater  than  0^. 

With  v,,  known,  the  Mach  number  in  region  n  is: 
un 


_2  -3 

1  +  1.3604  v  +  0.0962  v  -0.5127v 
1  -  0.6722  v  -  0.3278  v2 


(B-29) 


where 


-  T  (Vp I  -  i]  'B-30> 


The  ratio  of  static  pressure  to  total  pressure  in  region  n  is: 


-  1  +  T1  Mu. 


(B-31) 


and  the  pressure  coefficient  is 


(B-32) 


The  quantity  ptu  =  PtUl  with  Ptu^  given  by  Equation  (B-20)  or 

(B-26)  depending  on  whether  a  is  less  than  or  greater  than  0^. 
The  quantities  P00/Pt  and  qm/pt  are  given  by  Equations  (B-14) 
and  (B-15) . 

After  the  calculations  given  by  Equations  (B—  28 )  through 
(B-32)  have  been  repeated  sequentially  for  n  =  2,  3,...N,  the 
complete  pressure  distribution  on  the  upper  surface  of  the 
airfoil  has  been  determined. 


B. 2  CALCULATION  OF  CENTER-OF-PRESSURE  SHIFT  ON 
AN  AIRFOIL  DUE  TO  THICKNESS 


With  the  pressure  distribution  calculated  by  the  method 
described  in  the  preceding  section,  the  center-of-pressure 
location  on  the  airfoil  due  to  its  thickness  can  be  determined 
For  an  airfoil  with  zero  thickness,  the  center-of-pressure 
location  is  at  the  50-percent  chord.  Thus,  a  shift  due  to 
thickness  relative  to  this  point  will  be  calculated.  The 
center-of-pressure  is: 


(B-3 


where  Cm  is  the  sum  of  the  pitching  moments  produced  by  the 
various  regions  describing  the  airfoil,  taken  about  the  50- 
percent  chord,  and  CN  is  the  sum  of  the  normal  forces. 

Referring  to  Figure  B1  the  increment  in  normal-force 
coefficient  produced  by  the  n  region  is: 


>  -■ 


.•  .*  .•  -*  -»  -•  .*  .» 


where  l  is  some  nondimensionalizing  length.  Similarly,  the 

^  Vi 

pitching  moment  produced  by  the  n  region,  taken  about  the 
50-percent  chord,  c/2,  is 


ICm>n  -  4CP 


(xn  -  xn-l>  [l  <xn  +  xn-l>  1 
n  t  L2-  2c  J 


(x„  +  x_  .  ) 


(B-35) 


The  total  center-of-pressure  shift  due  to  the  thickness 
distribution  of  the  complete  airfoil  section  is  obtained  by 
summing  over  the  N  regions 


N 

I 

n=l 

N 

l  <CN). 
n=l 


(B-36) 
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APPENDIX  C 


EQUIVALENT  ANGLE  OF  ATTACK  FORMULATION 
FOR  LARGE  ANGLES  OF  ATTACK 

C.l  INTRODUCTION 

This  appendix  is  a  revised  and  expanded  version  of  Appendix 
C  of  reference  1.  The  purpose  of  this  appendix  is  to  derive  the 
equivalent  angle-of-attack  formulas  and  to  show  how  they  are  used 
(1)  to  develop  the  tables  to  be  incorporated  in  the  present  version 
of  PROGRAM  MISSILE  and  (2)  to  account  for  extrapolation  from  the 
measured  data  base.  In  particular,  it  is  useful  to  rederive  the 
formulas  in  order  to  handle  the  new  control  data  base  to  be 
obtained  in  the  future  which  was  not  available  at  the  time 
reference  1  was  written. 

The  equivalent  angle-of-attack  concept  is  based  on  the 
notion  that  the  aerodynamic  characteristics  of  an  all-movable  fin 
on  a  body  should  parallel  those  of  a  wing  alone*  with  the  same 
planform  and  the  same  Mach  number,  provided  that  the  effective 
incidence  of  the  flow  seen  by  the  fin  is  the  same  as  that  of  the 
wing  alone.  We  will  want  the  derivation  to  include  body  angle, 
deflection  and  span  effects  as  well  as  the  effects  of  vortices; 
i.e.,  we  want  to  find  an  <*eqi  f°r  i  such  that 

CNw^eq^  =  CNF^  (g)  tac»  4^  /  <5^  <5  2 , 5  3 , 6 4 ,  a/sm,  (Aaeq)Vil 

The  idea  for  aeg  was  suggested  by  the  earlier  slender-body 

work  of  Pitts,  Nielsen,  and  Kaattari  (ref.  43)  which  found  just 

such  behavior  for  the  linear  range.  However,  for  large  angles  of 

attack  a  nonlinear  definition  of  a  is  required.  Since  there  is 

eq 

no  unique  way  to  derive  a  nonlinear  formula  from  the  linear 


The  wing  alone  consists  of  two  opposing  fins  joined  at 
their  root  chords. 


result,  we  are  free  to  choose  our  approach  provided  that  it  gives 
the  linear  result  in  the  limit  of  small  angles.  Our  method  is 
based  on  the  use  of  average  velocity  components  seen  by  the  fin 
of  interest.  Those  velocity  components  are  put  together  to  give 
a  .  Extrapolation  from  the  measured  data  base  is  accomplished 
as  follows.  First,  the  equivalent  angle  of  attack,  oteg^, 
corresponding  to  the  vortex-free  value  of  fin  normal-force  coef¬ 
ficient,  CnF(Bj,  for  no  fin  deflection  is  obtained  from  the  wing- 
alone  normal-force  curve.  Then,  the  value  of  Aa^  corresponding 
to  changes  from  that  condition  is  determined  and  added  to  aGa  . 
The  new  value  of  a  is  used  together  with  the  wing-alone  curve 
to  obtain  the  corresponding  fin  normal-force  coefficient,  cNf(b)* 
The  procedure  is  illustrated  in  the  sketch  below.  The  details 
are  given  in  the  following  section. 


■F  (B) 


finish 


*F (B)  [start 


Illustration  of  equivalent  angle  of  attack  concept 


C.2  DERIVATION 


Consider  a  cruciform  wing-body  combination  with  the  x-axis 
rearward  along  the  body  axis,  the  y-axis  lateral  along  the  right 
horizontal  fin  (fin  4)  and  the  z-axis  vertical  along  the  upper 
fin  (fin  1) .  Let  the  combination  first  be  pitched  in  a  plane 
containing  the  free- stream  velocity  and  the  x-axis  by  an  angle 
ac.  Let  it  then  be  rolled  right  wing  down  by  angle  (Jig.  In  a 
plane  normal  to  the  body  axis  we  now  have  the  following  picture: 


V  sina 

00  C 


We  determine  first  the  flow  field  seen  by  the  fins  with  no 
fin  deflection  and  then  consider  the  effects  of  rotating  (deflec¬ 
ting)  the  fins  in  that  flow  field.  Referring  to  the  sketch  above 
we  see  that  the  components  of  the  free-stream  velocity  along  x,  y 


and  z  are  V00cosac,  -Voesinacsin^4  and  Vasina  ccos<J>4,  respectively. 
Hence,  the  component  of  free- stream  velocity  normal  to  the  plane 
of  fin  4,  Vn4,  is  V^sino^cos^ .  Thus,  the  angle  of  attack 
induced  on  fin  4  if  it  sees  the  free-stream  only  is 


tana 


'eq. 


=  tanaccos<}>4 


(C- 


Equation  (C-l)  does  not  yet  represent  the  actual  angle  of  attack 
induced  on  fin  4  with  no  fins  deflected  because  effects  of  the 
body,  sideslip  and  vortices  have  not  been  incorporated.  We  now 
consider  these  one  at  a  time. 


The  presence  of  the  body  affects  the  flow  in  the  crossflow 
plane  (wing-body  interference) .  We  account  for  this  effect  by 
multiplying  Vn4  by  the  factor  Kw,  i.e. 


Vn4 


free-stream+body 


=  KwV00sinaccos<J>4 


(C- 


For  favorable  interference,  Kw  >  1.  With  increasing  ac  and  M^, 

Kw  tends  to  decrease  (ref.  1) .  We  will  assume  that  Kw  is 
independent  of  <J>. 

Any  fin  for  which  <p  is  not  zero  is  sideslipping.  Spreiter 

and  Sacks  (ref.  45)  investigated  this  effect  and  found  that  the 

increment  in  the  fin  normal  force  was  proportional  to  the  product 

2 

of  a  and  8  for  the  fin,  i.e.  acsin$4cos<J>4 .  We  use  the  same  idea 
here  to  account  for  changes  in  Vn^  due  to  sideslip,  i.e. 


sideslip 


“  jk  %V„sin2acsin<f>4cos<fr4 


(C- 


Note  that  we  are  using  the  sine  definition  of  the  angle  of 
attack  (see  pg.  5  of  reference  44). 


C-4 


The  introduction  of  the  4/AR  coefficient  makes  the  slender  body 
value  of  independent  of  AR  for  delta  fins. 

To  account  for  the  effects  of  body  vortices  and  vortices 
generated  by  upstream  fins,  we  need  a  method  for  determining  the 
equivalent  angle  of  attack,  (Aaeq)v.,  induced  on  fin  i  by  the 
vortical  flow  field.  The  method  used  in  PROGRAM  MISSILE  is 
described  in  reference  1  and  is  given  next.  By  definition,  we 
have 


V-  =  tan(i“eq>v4 

vortices 


(04) 


The  method  assumes  that  the  vortices  are  rectilinear  and  parallel 

to  the  body  axis.  Hence,  (Vp4) vortices  ec3ua^  to  tlie  component 
of  the  free-stream  velocity  which  is  parallel  to  the  fin  chord 


and  we  can  write 


vortices 


=s  v  cosa  tan(Aa  ) 
w  c  eq  v . 


(C— 5) 


Using  relations  (02,3,5)  we  can  define  an  equivalent  angle 

/s 

of  attack,  <*eq4'  for  body,  sideslip,  and  vortex  effects  but  with¬ 
out  fin  deflection,  i.e. 


free- stream+bod 


sidesli 


vortices 


=  Kwtanaccos<|>4  K^tano^smo^sin^cos^  +  tan  (Aaeq) 

(06 

Note  that  the  sideslip  and  vortex  terms  in  equation  (06)  cannot 
increase  aeg^  beyond  ninety  degrees. 

To  account  for  the  change  in  equivalent  angle  of  attack  due 
to  fin  deflection,  we  define  a  new  quantity,  A such  that  the 


effect  on  fin  i  of  the  deflection  of  fin  j  is  given  by 


(Aa  )  . . 

eq  31 


-  Vi 


(C-7) 


where  is  to  be  obtained  experimentally  and  will  be  a  function 
of  6^.  Then  the  final  expression  for  the  equivalent  angle  of 
attack  of  fin  4  which  includes  fin  deflection  is 


a 


4 

BO  +  l 
#q4  j=l 


Vi 


(C-8) 


By  using  an  angle  addition  theorem  rather  than  a  tangent  addition 
theorem  we  allow  for  the  possibility  of  aeq4  exceeding  ninety 
degrees  due  to  large  fin  deflection.  In  the  definition  of  as 
incorporated  in  equations  (C-7)  and  (C-8) ,  we  have  assumed 
implicitly  that  any  vortex  effects  on  fin  i  are  not  changed  by 
deflecting  the  fins. 


We  can  generalize  equations  (C-6)  and  (C-8)  by  defining  <|>4 
to  be  the  bank  angle  of  fin  i  measured  positive  to  windward  from 
the  horizontal  plane.  The  horizontal  plane  is  defined  as  that 
plane  which  is  perpendicular  to  the  wind  plane  and  contains  the 
body  axis.  The  generalized  equations  are 


tana 


egi  *  Kwtanaccos*i  +  ^  K^tana^ina^in^cos^ 


and 


«<J. 


+  tan(ABeq) 


s +  jK 


(C-9) 


(C-10] 


Equations  (C-9)  and  (C-10)  represent  the  mathematical  formulation 
of  the  equivalent  angle  of  attack. 


LIST  OF  SYMBOLS 


A1'A2'A3 


CNw(B) 

CNC1 ,  CNC2 

CHMC1 , CHMC2 , 

CNTl ,CNT2 , .  . 

CHMT1,CHMT2, 


F1'F2'F3,F4 


body  radius  at  cylindrical  section 

empirical  coefficients  in  curve  fit  to  normal- 
force  curve;  see  Table  3 

aspect  ratio 

total  wing  span 

bending  moment  of  fin  about  root  chord,  posi¬ 
tive  tip  upward 

local  wing  chord 

local  wing  section  lift  coefficient 
wing  or  fin  root  chord 
wing  or  fin  tip  chord 
normal-force  coefficient,  N/qjSj^ 
pressure  coefficient,  (p  - 

root  bending-moment  coefficient,  BM/qooSR£r 

hinge-moment  coefficient,  HM/qooSRfr 

fin  normal-force  coefficient,  N/q^Sj^ 

fin  normal-force  coefficient,  N/q^Sj^ 

values  of  CNF  for  canard  fins  F^,F2#... 

.  values  of  CRM  for  canard  fins  F^,F2,... 

values  of  CNF  for  tail  fins  F1,F2,... 

.  values  of  CRM  for  tail  fins  F]_,  F2,... 

diameter  of  missile  body  on  cylindrical 
section;  2a 

fins  in  a  cruciform  arrangement  with  F^  at 
top;  numbering  is  clockwise  for  tail  fins  and 
counterclockwise  for  canard  fins 


Nv'.\ 


hinge  moment 

ratio  of  twice  fin-on-body  normal  force  at 
6  »  0  to  that  for  wing  alone  for  same  angle  of 
attack 

reference  length 

crossflow  Mach  number,  M  sina 

free- stream  Mach  number 

number  of  equally  spaced  fins  in  an  empennage; 
also  number  of  panels  on  chordwise  row  on  wing 
or  fin 

normal  force 

local  static  pressure 

free-stream  static  pressure 

free- stream  dynamic  pressure 

semi span  of  wing;  exposed  semi span  of  fin  on 
body 

semispan  from  body  axis  to  fin  tip  of  wing- 
body  combination 

reference  area 

wing-alone  planform  area 

maximum  thickness  of  a  wing 

thickness  ratio  of  airfoil  section 

coordinate  system  with  origin  at  leading  edge 
of  fin  root  chord  with  x  positive  rearward  and 
y  positive  laterally  towards  tip  looking  forward 
(see  Fig.  5) 

x,y  coordinates  of  fin  or  wing  center  of  pressure 
x,y  coordinates  of  fin  or  wing  center  of  pressure 
x  coordinate  of  wing  or  fin  leading  edge 


x  coordinate  of  fin  hinge  line 

axial  distance  of  orifice  in  front  of  fin 
trailing  edge 

axial  position  of  center  of  pressure  as  a 
fraction  of  fin  or  wing  root  chord 

value  of  x  for  fin  mounted  on  body  with  6=0 

value  of  x/cr  for  loading  on  fin  in  presence 
of  body  due  to  fin  deflection  6 

change  in  value  of  x  due  to  adding  thickness 
to  an  airfoil  or  wing 

value  of  x/cr  from  linear  theory  for  wing  with 
no  thickness 

value  of  x/cr  for  wing  at  a  =  45°  determined 
experimentally 

lateral  distance  of  body  vortex  from  body 
longitudinal  axis 

vertical  distance  of  body  vortex  above  x-y 
plane 

angle  of  attack 

angle  of  attack  of  missile  between  body  longi¬ 
tudinal  axis  and  free-stream  velocity;  different 
from  o  if  if  /  0 

equivalent  angle  of  attack  of  fin  with  no 
fin  deflections 

equivalent  angle  of  attack  of  fin  with  any 
or  all  fins  deflected 

angle  for  shock  detachment 

local  angle  of  attack  of  flow  approaching  fin 
leading  edge  in  plane  parallel  to  body  axis 
and  normal  to  fin  planform 


effective  aspect  ratio 
strength  of  body  vortex 

fin  deflection;  positive  trailing  edge  down 
for  horizontal  fins  and  trailing  edge  to  the 
right  for  vertical  fins 

fin  deflections  for  Fi*F2'F3»  and  F4'  resPectively 

fin  deflection  for  canard  fin 

fin  deflection  of  j***1  fin 

fin  deflection  for  tail  fin 

nose  semiapex  angle  of  wing  section  in  plane 
normal  to  wing  or  fin  planform  and  parallel  to 
root  chord 

value  of  0n  for  shock  detachment  in  two- 
dimensional  wedge  flow  at  a  *  0° 

wing  or  fin  taper  ratio,  ct/cr 

average  angle  of  attack  induced  on  fin  F.^  due 
to  unit  deflection  of  fin  Fj 

fraction  of  the  local  chord  of  the  fin  center- 
of-pressure  position  with  no  vortex  effects 

roll  angle  of  missile,  positive  clockwise  viewed 
from  rear 

roll  angle  of  fin,  measured  positive  clockwise 
from  positive  z  axis 

roll  angle  of  canard  fin  measured  positive 
clockwise  from  positive  y  axis 

roll  angle  of  tail  fin  measured  positive  clock¬ 
wise  from  positive  y  axis 
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